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SUMMARY 


INTRODUCTION 

The  present  work  is  an  improvement  and  application  of  previous  calcula- 
tional  methods  for  describing  building  wall  failure  and  movement  of  wall  and 
other  debris  by  airblast.  Both  wall  reactions  and  hydrodynamic  loading  are 
treated  by  means  of  lumped  parameters.  Classical  mechanics  using  springs 
friction  and  viscosity  is  applied  to  the  simulation  of  ground  interactions. 
Computerized  computations  of  even  large  building  collapse  and  subsequent  move¬ 
ment  of  debris  requires  considerably  less  time  than  would  be  needed  by  a 
finite  element  code  and  a  conventional  hydrocode. 

The  computation  system  is  presently  composed  of  two  manually  inter¬ 
facing  parts:  BRACOB  (Blast  Response  and  Collapse  of  Buildings)  which  cal¬ 
culates  loading  and  response  of  one  story  of  a  building  and  DEBRIS  which 
calculates  displacement  and  orientation  of  the  fragment  of  uniform  density. 
BRACOB  supplies  failure  time,  orientation  and  velocity  for  each  fragment,  pro¬ 
vided  the  user  defines  a  fragmentation  pattern.  The  blast  loading  on  the 
building  can  be  calculated  by  BRACOB  following  the  classical  prototype  or  can 
be  defined  arbitrarily  by  the  user.  The  blast  wind  in  DEBRIS  is  presently 
restricted  to  the  classical  form. 

THE  DEBRIS  MODEL  AND  OVER-THE-GROUND  TRANSPORT 

The  use  of  a  double-valued  spring  to  represent  the  ground  reaction  on 
the  fragment  corner  has  been  extended  to  both  the  vertical  and  horizontal 
directions.  A  double-valued  spring  constant  results  in  the  dissipation  of 
energy  apparently  more  effectively  than  simple  Coulomb  friction  and  velocity- 
proportional  viscosity.  It  was  found  friction  and/or  viscosity  without  horizon¬ 
tally-acting  ground  springs  led  to  high  speed  rotation  of  the  fragment  which  in 
turn  greatly  reduced  the  relative  speed  over  the  ground  of  each  corner  and  thus 
also  reduced  the  viscous  loss.  When  the  frictional  or  viscous  forces  were  in¬ 
creased  still  further  in  an  attempt  to  compensate,  bouncing  and  loss  of  the 
contact  with  the  ground  occurred.  In  any  case,  when  only  friction  and  viscous 
losses  are  used,  the  DEBRIS  code  calculates  transport  of  some  of  the  light¬ 
weight  concrete  cubes  that  is  excessive  when  compared  to  observations  at  the 
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DIAL  PACK  event  (23  July  1970).*  This  lack  of  agreement  between  simulations 
and  observations  occurs  at  incident  peak  overpressures  of  50  psi  and  greater. 
However,  reducing  the  stiffness  of  the  vertical  ground  spring  and  adding  a 
double-valued  horizontally  acting  ground  spring  has  made  it  possible  to  simulate 
the  total  distance  transported  by  1  kt  blasts  over  the  range  of  incident  pres¬ 
sures  from  30  to  100  psi,  as  can  be  seen  in  Figure  S-l.  (In  the  Figure,  K^ 
and  refer  to  the  two  values  of  the  spring  constant  for  the  vertically  acting 
ground  spring  and  K  ^  and  K  ^  r,?Present  values  for  the  horizontally  acting 
ground  spring.  The  values  u  and  n  stand  for  the  Coulomb  friction  and  viscosity, 
respectively).  Moreover,  the  simulated  motion  is  a  tumbling  in  which  the  cube 
never  rises  more  than  1.5  inches  above  the  ground.  The  height  of  rise  during 
tumbling  is  controlled  by  the  value  of  the  vertical  ground  spring.  Increasing 
the  stiffness  by  a  factor  of  ten  can  increase  the  rise  by  an  order  of  magnitude. 
The  distance  transported  is  also  increased  substantially. 


Figure  S-2  shows  the  gradual  braking  action  of  the  DIAL  PACK  simulation 
containing  both  the  vertical  and  horizontal  ground  springs.  The  two  ordinates 
show  horizontal  speed  and  angular  speed  of  the  tumbling  cube  originally  placed 
at  the  50  psi  contour. 

Over  the  range  of  incident  over  pressure  30  to  100  psi  drag  or  dynamic 
pressure  forces  in  the  blast  completely  dominate  diffraction  forces  in  the 
transport  of  the  DIAL  PACK  cube.  At  15  psi  drag  is  much  weaker  than  at  30  psi 
and  above  so  that  an  estimate  of  momentum  gained  through  shock  diffraction  must 
be  added  to  the  simulated  cube  before  agreement  with  observation  is  obtained. 


The  quantity  At,  which  is  identified  in  Figure  S-2,  represents  the  size 
of  the  time  mesh  and  has  been  found  to  be  important  although  no  quantitative 
relation  has  been  derived,  as  yet. 

The  values  of  the  simulation  parameters  entered  into  Figure  S-l  are  close 
to  the  optimum  for  simulation  of  the  transport  of  the  DIAL  PACK  cube  which  was 
initially  resting  on  the  ground  surface.  Other  experiments  in  which  objects  are 
dropped  from  moving  trucks  at  a  height  of  approximately  2.6  feet  above  the 
ground  can  not  be  simulated  with  these  "optimum"  values.  In  particular,  the 
horizontal  spring  must  be  weakened  by  two  orders  of  magnitude  in  order  to  achieve 


*  Concrete  cubes  weighing  65  lbs.  and  measuring  1  foot  on  edge  were  placed  on 
the  ground  at  various  distances  from  Ground  Zero  and  total  distance  moved  during 
the  explosion  recorded.  Explosion  yield  was  approximately  500  tons  of  TNT. 
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Figure  S-2 .  Braking  of  Concrete  Cube  at  DIAL  PACK. 


realistic  simulations  of  distances  transported  in  these  experiments  from  moving 
trucks.  When  this  change  is  made  however,  some  of  the  observations  can  be 
simulated,  specifically,  the  distance  travelled  by  the  dropped  object  is  in¬ 
dependent  of  object  density  and  is  within  one  standard  error  of  estimate  of 
observed  values.  Whether  or  not  the  need  to  change  the  value  of  the  horizontal 
spring  is  related  to  the  differences  in  the  surfaces  over  which  the  two  kinds 
of  experiments  were  conducted  is  not  known.  The  photographs  from  DIAL  PACK  show 
a  rough  surface  covered  with  tufts  of  grass  while  the  experiments  from  moving 
truer s  were  conducted  over  a  "graded  airstrip  built  on  an  alluvial  plain". 

Using  the  optimum  parameter  values,  as  described  above,  distance  the 
simulated  DIAL  PACK  cubes  are  transported  almost  scales  with  the  cube  root  of 
yield. 


COLLAPSE  OF  THE  LANDIS  HOSPITAL 


The  exterior  walls  of  the  Henry  R.  Landis  State  Hospital,  Philadelphia, 
one  of  the  buildings  in  the  National  Shelter  Survey,  consist  of  two  tiers  of 
masonry  each  capable  of  arching  in  a  nearly  rigid  frame.  A  cross  sectional 
view  appears  in  Figure  S-3.  The  outer  tier  of  clay  brick  arches  (one-way) 
between  steel  supports;  the  inner  tier  of  concrete  masonry  units  arches  (two- 
way)  in  a  reinforced  concrete  frame.  In  the  analysis  both  frames  are  treated 
as  rigid.  The  two  tiers  respond  to  air  blast  in  parallel  but  because  of  the 
extreme  thinness  of  each  tier  the  incipient  collapse  overpressure  of  the  ex¬ 
terior  walls  in  less  than  1  psi.*  The  differences  between  tne  simulated  be¬ 
haviors  under  1  and  30  psi  blast  waves  are  entirely  in  the  debris  histories. 

Peak  Free-Field  Overpressure  =  1  psi 

Interior  wall  behavior  has  a  small  but  measureable  effect  on  exterior 
wall  response.  Under  1  psi  blast  wave  loading  the  pressure  of  interior  par¬ 
titions  may  delay  failure  of  the  front  exterior  walls  by  as  much  as  10  ms  and 
may  reduce  the  departure  speed  of  the  fragments  from  approximately  11  to  6  fps. 
Figure  S-4  is  the  floor  plan  that  was  analyzed  in  this  study.  The  two  coordinates 
(marked  I  and  J  in  the  Figure)  locate  corners  and  openings  for  the  computer 
calculation.  It  is  estimated  that  the  front  or  south  elevation  will  fail  under 
1  psi  blast  50  to  60  ms  after  blast  arrival  and  initial  center  of  mass  fragment 
speed  will  be  in  the  range  6.2  to  7.8  fps.  All  interior  east-west  walls  will  be 
breached  almost  upon  blast  arrival,  except  the  walls  at  grid  line  I  =  11,  where 
the  build-up  of  pressure  against  the  surface  facing  the  origin  of  the  blast 
must  await  filling  of  a  large  volume.  However  pressure  buildups  on  the  opposite 
(north)  face  is  also  slow  because  the  windows  in  the  north  elevation  of  the 
building  lie  in  the  wake  of  the  blast.  It  is  estimated  the  interior  partition 
line  at  I  =  11  fails  sometime  between  50  and  70  ms  after  blast  arrival.  This 
failure  is  the  prelude  to  outward  failure  of  the  north  exterior  walls  at  170  ms 
with  a  departure  speed  of  approximately  3.4  fps.  At  this  point  the  blast  from 
a  1  Mt  weapon  still  has  several  seconds  of  life  left  and  further  downwind  trans- 

*  Collapse  overpressure  is  stated  in  terms  of  peak  free-field  overpressure  in 
a  classical  blast  wave  oriented  head-on  to  one  facade  of  the  building. 
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port  of  debris  occurs. 

At  1  psi  failure  of  the  side  (east  and  west)  walls  is  marginal  and  results 
of  calculations  depend  critically  on  assumptions  of  behavior  of  interior  par¬ 
titions.  Most  probably  sidewalls  fail  inwardly  at  very  low  speeds. 

Peak  Free-Field  Overpressure  =  30  psi 

Under  30  psi  blast  loading  the  south  exterior  walls  fail  at  approximately 
6  ms  with  departure  speeds  in  the  neighborhood  of  110  fps.  Simple  room-filling 
through  the  south  windows  quickly  breaches  the  east-west  interior  partitions  in 
the  front  tier  of  rooms  and  the  prompt  failure  of  the  south  exterior  walls  in¬ 
sures  the  presence  of  loading  in  excess  of  10  psi  on  the  south  face  of  interior 
partition  line  I  =  11  well  before  40  ms  when  the  loading  of  the  north  exterior 
facade  begins.  Blast  filling  through  the  north  windows  then  has  no  chance  of  de¬ 
laying  failure  of  the  interior  wall  line  at  1  =  11.  Consequently  the  north 
exterior  walls  fail  outward  at  approximately  62  ms  with  departure  speeds  of  about 
8  fps. 

Side  walls  at  30  psi  will  most  likely  fail  inward.  For  example, 
southernmost  sidewalls  fail  at  23  and  17  ms  with  speeds  of  27  and  40  fps,  res¬ 
pectively.  Due  to  the  flexure  of  these  walls,  however,  much  of  their  debris  will 
be  deposited  outside  the  building. 

FRAGMENT  DISTRIBUTION,  LANDIS  HOSPITAL 

Even  using  loss  parameters  found  to  be  appropriate  for  a  smooth  hard 
surface  the  further  translation  of  collapsed  wall  fragments  by  the  1  psi  blast 
wind  does  not  appear  significant.  Fragments  from  the  south  exterior  wall  would 
appear  to  be  distributed  over  the  floor  in  an  area  reaching  from  2  to  9  ft.  in¬ 
side  the  initial  wall  line.  Height  of  bounce  is  never  more  than  a  few  inches, 
a  circumstance  related  to  the  relatively  slow  rotation  of  the  fragments  upon 
their  departure  from  the  wall.  Fragments  from  the  north  exterior  wall  move 
even  less  from  their  points  of  impact  than  does  the  south  wall  debris,  but  since 
these  fragments  descend  from  positions  in  a  six  story  wall  to  the  ground,  the 
total  horizontal  transport  ranges  from  3  to  15  ft. 

The  fragment  disposition  at  30  psi  is  dramatically  different.  The 
first  fragments  from  the  south  wall  reach  the  floor  approximately  41  ms  after 
blast  arrival  at  a  distance  of  3  ft.  from  the  wall  line.  Aejuming  these  frag¬ 
ments  remain  intact,  they  then  bounce  to  a  height  of  approximately  2.5  ft. 
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above  the  floor  and  travel  72  ft.  downstream.  Such  behavior  places  south  wall 
debris  at  the  north  exterior  wall  line  at  370  ms.  During  this  transport 
collision  with  interior  wall  debris  will  be  highly  unlikely  since  such  light 
material  as  is  found  in  the  interior  partitions  will  surely  be  swept  ahead  of  the 
masonry  material. 

After  their  first  impact  with  the  concrete  floor  north  wall  fragments 
<  may  be  further  fragmented.  Treating  this  debris  as  individual  masonry  units 

moving  initially  with  the  speed  of  the  center  of  fragment  mass  the  horizontal 
transport  is  essentially  the  same  as  described  in  the  preceding  paragraph  ex¬ 
cept  the  motion  is  a  tumbling.  Under  this  assumption  south  wall  debris  reaches 
the  north  facade  at  340  ms,  long  after  the  north  wall  fragments  have  departed. 

Despite  their  high  initial  speed  and  early  departure,  the  south  wall 
fragments  will  not  necessarily  overtake  the  fragments  from  the  north  wall. 
Trajectories  differ  markedly.  By  170  ms  a  near  classical  blast  wave  is 
blowing  through  the  building.  The  rapidly  rotating  south  wall  fragments  are 
moved  upward  or  sidewise  relatively  little.  Rear  wall  fragments,  on  the  other 
hand,  keep  for  a  much  longer  time  an  orientation  that  permits  lofting  of  top 
wall  fragments  or  side  motion  in  side  fragments.  The  consequence  is  a  far 
longer  trajectory  in  rear  wall  than  in  front  wall  fragments,  as  is  illustrated 
in  Figure  S-5.  The  small  arrows  in  the  Figure  show  the  attitudes  of  the  frag¬ 
ments.  In  neither  case  do  the  fragments  travel  far  over  the  ground.  Fragments 
from  the  first  (ground)  floor  of  the  south  wall  strike  the  ground  in  the  wake 
of  the  building  while  blast  winds  are  still  blowing,  yet  travel  only  55  ft. 
after  touching  down.  Comparable  fragments  from  the  sixth  floor  move  only  10 
ft.  after  striking  the  gr'-rtv'. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  techniques  developed  and  elucidated  in  this  research  are  certainly 
capable  of  describing  gross  features  of  debris  distribution  from  some  collap¬ 
sing  building  walls;  e.g.,  in  the  Landis  Hospital  the  dramatic  differences  be¬ 
tween  the  effects  of  1  and  30  psi  blast  waves,  the  relative  timing  of  collapse 
events  and  debris  movements.  But  some  finer  points  await  improvement  in  the 
tools  themselves;  e.g.,  the  conditions  for  removal  of  interior  partitions  must 
be  made  more  precise  and  incorporated  into  the  automatic  calculations,  data  to 


10 


.20  x  10  dyne/cm 


o 

o 


Figure  S-5.  Lofting  of  Ground  Floor  Top  Wall  Fragments,  Landis  Hospital,  p  =30  psi. 


describe  debris  transport  over  a  concrete  surface  must  be  found  before  full 
confidence  can  be  placed  in  some  of  the  present  calculations,  and  finally  it 
must  be  pointed  out  that  the  case  treated  in  this  research,  i.e.,  a  weak  walled 
building  with  a  relatively  strong  frame  presents  a  relatively  simple  debris 
distribution  problem.  Weak  frame  buildings,  tilt-up  structures,  and  certain 
massive  structures  supported  by  load-bearing  masonry  walls  contain  complications 
that  may  perturb  the  debris  distribution  described  here. 
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DEBRIS  DISTRIBUTION  AS  A  PARAMETER 
IN  BLAST/FIRE  INTERACTION 


PART  I 


INTRODUCTION 

Previous  civil-defense-inspired  work  contains  estimations  of  the 
quantity  and  distribution  of  structural  debris  stemming  from  the  effects 
of  air  blast  on  buildings . 1 » 2 »  ■ » 4 » 5  The  estimates  were  apparently  made 
to  determine  accessibility  and  the  need  for  debris  clearance.  The  present 
research,  done  in  support  of  studies  on  the  blast/fire  interaction,  has 
been  designed  to  provide  more  detailed  descriptions  of  temporal  and  spa¬ 
tial  debris  configurations  than  have  previous  civil-defense  studies.  Pre¬ 
viously,  rather  detailed  efforts  have  determined  debris  trajectories  in 
connection  with  U.S.  missile  defense. 6 * 7 > 8  Some  of  the  missile  defense 
work  has  been  most  useful  to  the  present  research,  other  parts  that  deal 
with  small,  unrealistically  constructed  models  have  not  been  used. 

The  present  effort  by  no  means  provides  all  the  answers  required  to 
completely  understand  the  effects  of  debris  on  the  interaction  of  blast 
and  fire.  The  focus  is  on  wall  debris,  excluding  room  contents  such  as 
furniture.  The  emphasis  is  on  exterior  walls.  Other  facets  of  the  blast/ 
fire  problem  are  important  and  can  be  examined. 

The  method  adopted  here  is  based  on  two  simulation  mod  Is  embodied 
in  the  computer  codes  BRACOB  and  DEBRIS.  BRACOB  is  a  system  of  codes 
that  compute  the  simultaneous  response  of  all  the  walls  on  one  floor  of 
a  building.  At  each  step,  it  determines  the  net  loading  on  each  wall, 
finds  the  motion  of  the  wall,  and  adjusts  the  loading  when  wall  failures 
change  the  geometry  of  the  building.  The  initial  loading  may  arise  in  a 
classical  blast  wave  or  in  an  arbitrarily  defined  pulse,  such  as  is  used 
in  taking  into  account  the  blast  interference  by  other  buildings. 

DEBRIS  determines  trajectories  of  single  fragments  through  the  air 
or  along  the  ground  in  response  to  classical  blast  waves.  Not  treated 


in  DEBRIS  is  the  impact  of  one  fragment  on  another  and  the  effects  of 
fragment  breakup.  Although  BRACOB  supplies  time  of  wa1 1  failure*  as 
well  as  wall  speed  and  displacement  at  the  moment  of  failure  ,  the  pat¬ 
tern  of  wall  fragmentation  must  be  defined  by  the  system  user.  For¬ 
tunately,  a  substantial  body  of  observation  exists--particularly  from 
atmospiieric  nuclear  weapons  tests0  and  from  tunnel  b  lasts^  n  >  •  -1  »  '  - — to 
help  define  crack  patterns. 

The  single-degree-of-freedom  engineering  analysis  of  wall  response 
used  in  BRACOB  is  from  Wiehie  and  Bockholt 1 0 » 1 4  and  has  been  applied  suc¬ 
cessfully  to  civil  defense  prob l-mis .  1  5 > 1  6  > 1  7 > 1  8 * 1  °  The  computer  simula¬ 
tion,  DEBRIS,  is  drawn  from  engineering  methods  of  aerodynamics  and  from 
classical  mechanics.20  One  of  the  goals  of  the  present  research  was  the 
continued  calibration  of  the  model  against  the  limited  data  available. 

This  report  falls  into  two  parts:  (1)  an  exposition  of  further  im¬ 
provements  in  the  computational  simulations  and  (2)  simulations  of  col¬ 
lapse  and  debris  transport  within  an  actual  building. 

OVER-TIIE-GROUND  TRANSPORT 

Any  simulation  of  dtoris  distribution  by  megaton  nuclear  blast  must 
accurately  estimate  the  interaction  of  a  moving  fragment  with  the  ground 
or  other  large  surfaces.  Except  in  buildings  over  three  stories  in 
height,  wall  debris  will  reach  the  ground  relatively  early  in  the  positive 
dynamic  pressure  phase.  Even  in  much  taller  buildings,  wall  debris  will 
reach  the  ground  with  a  very  high  horizontal  velocity  component.  How  the 
ground  brakes  the  blast-driven  fragment  is  crucial  to  determining  the 
spread  of  fragments. 

Previous  Work 

Fortunately,  observations  of  blast-driven  objects  exist.  At  opera¬ 
tion  DIAL  PACK  in  July  1970,  severa’  < 1 ight-weight)  concrete-filled  ply¬ 
wood  cubes  were  placed  on  the  ground  at  various  ranges  from  ground  zero 

*"Failure"  is  defined  as  the  moment  when  the  component  ceases  to  fulfill 
its  structural  function,  e.g.,  support  a  vertical  load.  "Collapse"  will  be 
used  to  mean  the  less  well  defined  point  at  which  the  component  no  longer 
obstructs  air  flow. 
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and  driven  downrange  by  the  blast  from  the  detonation  of  approximately 
500  tons  of  high  explosive  placed  on  the  surface.  A  postshot  photograph 
of  these  cubes  is  reproduced  in  Figure  1.  Incident  peak  overpressures 
at  the  cubes  varied  from  15  to  100  psi  (103  to  689  kPa) .  The  cubes  were 
one  foot  on  an  edge21  and  only  their  displacements  by  the  blast  were  re¬ 
ported;  no  information  on  trajectories  has  been  found,  but  distances 
travelled,  weights  of  cubes,  and  initial  peak  overpressures  are  shown 
in  Table  1,® 

Fletcl  r  and  Bowen22’23  have  sh jwn  experimentally  that  the  distance 
required  to  brake  stones  and  concrete  blocks  to  a  stop  depends  only  on 
the  initial  speed.  In  these  tests,  the  objects  were  thrown  out  of  a 
moving  truck  at  a  height  of  approximately  2.6  ft  above  a  graded  airstrip 
built  on  an  alluvial  plain;  results  for  stones  and  concrete  masonry  blocks 
were  reported  in  two  recursion  relations: 

log^Q  s  =  -1.4664  +  1.9004  log^  V 
log  t  =  -1.1420  +  0.9004  log  1Q  V 

where  s  =  distance  (ft)  from  point  of  release  to  rest  position,  t  =  time 
(s)  from  time  of  release  to  stopping  time,  and  V  =  horizontal  speed  (ft/s) 
of  truck  at  moment  of  release.  These  regression  relations  are  nearly 
equivalent  mathematically  to  the  existence  of  a  constant  force,  F,  pro¬ 
portional  to  the  mass,  m,  and  to  a  fractional  power  of  the  initial  speed, 
V: 


F  = 


14.63  V 


.0996 


m 


acting  through  the  duration  of  the  motion  of  the  object.  The  stones 
weighed  from  1.3  to  116  lb;  tfu  masonry  blocks  were  both  hollow  (25  to 
32.5  lb)  and  concrete-filled  (55  to  56.5  lb).  Standard  error  in  the  es¬ 
timate  is  +41.2  percent.  The  attitude  of  the  object  when  dropped  was 
randomly  varied,  and,  no  doubt,  this  randomness  contributed  to  the  large 
error  of  estimate. 
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Table  1  .  Weight  and  Transport  of  Tumbling  Blocks  (from  Ref.  8). 


Tumbling 

Block 

Number 

Tumbling 

Block  Initial 
Overpressure 
Position 
(psi) 

Actual 

Weight 

(Pounds) 

Approximate 

Radial 

Transport 

Distance 

(ft) 

1 

15 

64.5 

6.7 

2 

15 

65.5 

7.3 

>  mean  7.5 

3 

15 

64.0 

8.5 

/ 

4 

30 

63.0 

36  ' 

I 

5 

30 

64.0 

39 

>  mean  39 . 3 

6 

30 

61.0 

43  J 

7 

50 

63.5 

103  ' 

8 

50 

63.5 

150 

>  mean  146.0 

9 

50 

61.5 

184 

/ 

10  100  Destroyed 

11  100  63.5  230  >  mean  273.0 

12  100  61.5  317 
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Using  aeronautical  engineering  measurements  of  wind  drag  and  nodeling 
ground  forces  by  double-valued,  vertically  acting  springs  and  coulomb 
friction,  the  computer  code  DEBRIS  has  successfully  simulated  the  data  in 
Table  1  over  partial  pressure  ranges  in  previous  work.90  A  major  task  of 
the  present  effort  was  to  devise  satisfactory  simulations  for  the  whole 
range  of  data  stemming  from  the  tumbling  cubes.  Techniques  that  produce 
excellent  results  at  low  incident  overpressures  fail  at  higher  pressures. 
Because  trajectory  or  displacement  versus  time  information  has  not  been 
available,  we  do  not  know  whether  the  shortcomings  of  the  previous  model 
lie  in  the  acceleration  or  deceleration  phases  of  the  transport.  How¬ 
ever,  while  wind  drag  is  a  subject  for  which  a  great  deal  of  literature 
exists,  relative  paucity  characterizes  the  engineering  literature  of  tum¬ 
bling  deceleration.  The  next  section  contains  the  results  of  the  research 
to  find  an  engineering  simulation  of  ground-fragment  interaction. 


Improvements  in  the  Model  DEBRIS 

With  a  friction  force  opposing  horizontal  motion  and  a  two-valued 
spring  exerting  vertical  ground  resistance,  simulated  air-blast  transport 
of  concrete  objects  generally  agrees  with  observations  at  peak  free-field 
overpressures  of  15  psi  and  below.  The  simulated  observations  were  those 
of  the  tumbling  blocks  at  15  psi  ai  DIAL  PACK8  and  the  collapse  of  a 
concrete  masonry  building  situated  at  the  9-psi  contour  at  PRAIRIE  FLAT. 1 
The  frictional  force,  F,  is  proportional  to  a  coefficient  of  friction,  p, 
and  the  force  between  the  object  and  the  ground,  i.e., 


Ff  =  ^Fz 


where  z  is  the  coordinate  in  the  direction  perpendicular  to  the  ground 

and  F^  is  the  force  in  the  ground  reaction  spring.  Force,  F ,  is  exerted 

in  the  opposite  direction  to  the  horizontal  motion.  As  the  corner  of  the 

£ 

object  moves  below  the  ground  surface  under  the  influence  of  gravity, 


If  the  ground  and  fragment  surface  were  in  intimate  contact,  air  blast 
static  pressure  would  furnish  a  second  downwind  force  against  the  verti¬ 
cal  spring.  However,  a  reai  ground  surface  will  probably  admit  static  air 
blast  pressure  underneath  the  fragment.  In  this  case,  the  net  downward 
air  blast  pressure  is  zero. 
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the  ground  reaction  spring  exerts  an  upward  force: 


F  =  -k  , z ,  when  z  <  0 
z  gl 


When  the  corner  is  still  below  the  surface  but  moving  upward, 


F 

z 


-k  „z 
g2 


so  that,  if  Az  is  the  maximum  downward  penetration,  the  energy  loss  in  a 
single  impact  is 


AE  =  i(k  ,  -  k  0)Az2 
2  gl  g2 


Handbook  values  of  coefficient  of  friction  vary  from  0.02  to  4.0, 
but  the  majority  of  values  lie  in  the  range  0  to  0.5. 24  None  of  the  ma¬ 
terials  listed  in  the  handbook  resemble  masonry  on  soil  in  gross  physical 
properties.  In  simulations  of  the  DIAL  PACK  tumbling  cube,  a  value  of 
p  =  0.4  predicts  nearly  the  reported  transport  at  15  psi,  but  at  50  and 
100  psi,  the  simulated  transport  distances  greatly  exceed  those  observed. 
From  Table  2,  increasing  the  coefficient  of  friction  to  0.5  improves  the 
simulation  at  the  high  original  overpressures,  but  results  in  less  than 
observed  transport  at  low  pressures.  Furthermore,  when  u  =  0.5  and  when 
original  pressure  is  100  psi,  the  transport  is  far  larger  than  observed; 
in  fact,  at  Pgo  =  100  psi,  the  coefficient  of  friction  should  be  in  the 
neighborhood  of  0.75,  if  reproduction  of  observed  transport  is  to  be 
achieved . 

There  is  another  unrealistic  feature  of  the  motion  simulated  by  means 
of  coulomb  friction  with  coefficients  in  the  range  0  to  1.  In  Figure  2, 
the  angle  of  tilt  of  the  cube  during  simulated  transport  from  an  initial 
position  at  Pso  =  100  psi  is  plotted  as  a  function  of  downrange  distance, 
y.  The  simulated  cube  never  loses  contact  with  the  surface  over  which 
it  moves  and  never  tilts  significantly.  Such  movement  can  hardly  be  de¬ 
scribed  as  "tumbling."  Even  from  the  photograph  in  Figure  1,  it  is 
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Dynamic  Pressure 


Figure  2.  Angle  of  Tilt  of  Concrete  Cube  at  DIAL  PACK  (friction  only). 


apparent  there  could  have  been  no  purely  s]iding  motion  of  the  concrete 
cubes  over  the  ground  at  DIAL  PACK. 

Simulated  tumbling  can  be  achieved  by  increasing  the  value  of  p. 

The  high  frictional  force  on  the  bottom  of  the  initially  sliding  cube 
causes  it  to  rotate,  followed  by  compression  of  the  ground  spring  and 
ending  with  upward  motion  of  the  cube  and  loss  of  contact  between  ground 
and  cube.  The  subsequent  bouncing  motion  is  capable  of  permitting  a  much 
greater  transport  distance  than  is  achieved  by  a  purely  sliding  contact. 
Under  these  circumstances,  the  downrange  distance  moved  depends  on  the 
value  of  the  vertical  spring  constant,  which  will  be  evident  later.  Even 
setting  the  restoring  spring  constant,  kg2,  equal  to  zero  does  not  prevent 
upward  motion  because  rotation  may  cause  the  corners  ot  the  cube  to  move 
into  the  ground  at  the  same  time  the  center  of  mass  is  moving  upward, 
and  the  "braking"  spring,  may  contribute  to  upward  acceleration. 

Some  results  of  high  value  for  the  coefficient  of  friction  are  shown  in 
Table  2  for  an  initial  location  at  the  30  psi  contour.  The  asterisk  in 
the  table  indicates  that  the  corresponding  motion  includes  airborne  seg¬ 
ments  or  tumbling.  Thus,  increasing  p  from  1.0  to  1.2  shortens  the  dis¬ 
tance  transported,  but  increasing  p  still  further  allows  for  greater 
transport  than  when  p  is  small.  Some  other  characteristics  of  the  motion 
of  this  simulation  at  p  =  30  psi  are  represented  in  Figures  3  and  4, 

i.e.,  height  of  center  of  mass  above  ground  and  angle  of  tilt,  both  as 
functions  of  the  downrange  distance  transported.  Although  the  total 
transport  in  the  case  of  p  =  1.4  and  p  =  30  psi  appears  nearly  correct, 

the  same  value  of  p  in  the  case  of  p  =  100  psi  yields  a  distance  trans- 

so 

ported  that  is  much  larger  than  observed.  Hence,  even  though  tumbling 
motion  may  be  simulated  with  friction  only,  the  coefficient  of  friction 
must  still  be  varied  with  incident  overpressure. 

Viscosity .  A  viscous  force,  i.e.,  resistance  proportional  to 
speed  or  some  power  of  speed,  seems  to  provide  braking  that  increases 
with  peak  overpressure.  It  is  easy  to  estimate  a  value  of  a  coefficient 
of  viscosity,  q- 
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Figure  4.  Simulated  Rotation  of  Concrete  Cube  at  DIAL  PACK  (friction  only) 


f  =  ny 


where  F  is  the  force  resisting  horizontal  motion  and  y  is  the  time  deriva¬ 
tive  of  the  horizontal  component  of  the  object's  motion. 

A  zoiuth  order  numerical  approximation  to  p  can  be  found  by  treating 
the  viscous  effect  as  a  correction- 

'«/  ft dt 

in  which  yc  represents  the  observed  displacement  of  the  center  of  mass, 
acceleration,  a,  is  calculated  from  wind  and  friction  forces  only,  yD 
is  intial  fragment  speed  (if  any),  and  M  =  cube  mass.  At  DIAL  PACK  y0  = 

0,  mass  M  =  65  lb  and  the  second  term  on  the  right  can  be  estimated  from 
Table  2.  The  equation  can  be  solved  as  follows: 


'A  111  +/  / 


a  dt  dt 


y  (calculated)  -  y  (observed ) 
/yd 


M 


Table  2  shows  for  overpressure  of  100  psi  that  the  calculated 
value  of  y  is  503  ft  when  p  =  0.5  and  no  viscosity  is  used.  The  dif¬ 
ference,  therefore,  between  calculated  and  observed  distances  on  the 
right  side  of  the  equation  above  is  estimated  to  be  230  ft;  and  the  dura¬ 
tion  of  simulated  motion  is  the  integral  representing  the  area  under  the 
calculated  displacement  curve  is  approximately  1901  ft-s.  Because  the 
mass  of  the  cube  is  reportedly  65  lbs,  the  first  estimate  of  n  becomes 
7.86  lbs/s  or  3567  g/s. 

When  the  displacement  history  is  recalculated  using  n  =  3000  g/s, 
the  plot  of  the  trajectory  in  a  plane  perpendicular  to  the  ground  (i.e., 
in  the  y-z  plane)  is  presented  in  Figure  5.  Total  displacement  is  overly 
large,  but  the  reason  is  clear  from  the  illustration:  the  cube  bounces 
so  much  that  it  is  seldom  in  contact  with  the  ground  where  the  viscous 
and  frictional  forces  operate. 
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•35  x  10  dyne/cm 


Simulated  Traiectorv  of  Conrrpt-P 


Another  revealing  feature  of  the  simulation  with  q  =  3000  g/s  is 
demonstrated  in  Figure  6,  a  plot  of  tilt  angle  against  downrange  distance, 
yc .  Clearly,  the  cube  is  spinning  at  a  high  rate  (~  18  rps) .  In  tact, 
detailed  examination  of  the  computer  output  shows  that  the  over-the-ground 
speed  of  the  lower  corners  is  much  less  than  the  downrange  speed  of  the 
center  of  mass.  Inclusion  of  the  viscous  force  appears  to  be  self- 
defeating:  the  net  result  has  been  to  reduce  the  braking  force. 

The  general  validity  of  this  observation  may  be  verified  by  computing 
a  trajectory  without  permitting  a  torque  about  the  internal  a-axis.  Fig¬ 
ure  7  is  two  such  displacement  histories  corresponding  to  two  different 
values  of  viscosity  for  an  original  location  on  the  100-psi  contour. 
(Originally  the  internal  a-axis  coincides  with  the  exlernal  x-axis;  the 
blast  direction  is  in  the  positive  y-direction,  so  that  rotation  caused 
by  contact  with  the  ground  is  about  the  a-axis. ?0)  Clearly,  the  succes¬ 
sive  approximations  outlined  above  will  converge  on  a  successful  value 
of  n  as  long  as  no  torque  caused  by  ground  forces  are  allowed.  Such  a 
restriction  is  unrealistic. 

Although  it  is  meaningless  to  continue  the  series  of  successive  ap¬ 
proximation  using  the  results  of  the  simulation  with  n  =  3000  g/s,  trial 
and  error  using  several  values  of  n ,  but  without  the  restriction  on 
torque,  produce  the  results  tabulated  in  Table  3.  In  all  these  cases 
coefficients  of  friction,  p,  is  0.5,  but  three  different  values 
of  q  are  represented.  With  two  exceptions,  all  simulated  motions  slide, 
including  usually  slight  rocking  back  and  forth,  as  illustrated  by  a  plot 
of  tilt  versus  ground  range  in  Figure  8.  When  PSQ  =  100  psi  and  p  = 

1500  g/s,  transport  is  not  by  sliding;  rather  the  simulated  cube  bounces. 
The  height  of  the  center  of  mass  is  drawn  in  Figure  9  as  a  function  of 
ground  range.  In  this  case,  as  well  as  the  case  shown  in  Figure  5,  the 
cube  is  spinning  (reaching  a  maximum  rate  of  approximately  3  rps),  but 
the  rate  of  spin  of  the  cube  represented  by  Figure  9  is  not  high  enough  to 
reduce  the  viscous  force  to  a  negligible  quantity. 

Displacements  reported  in  Table  3  suggest  that  the  viscous  force 
appropriate  to  this  simulation  may  not  be  proportional  to  the  first  power 
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Simulated  Tilt  of  Concrete  Cube  at  DIAL  PACK  (high  viscosity) 


Dynamic  Pressure 
Extinction 


of  Ground-Induced  Torque  Set  to 


Pso^ 

(psi) 

n 

15 

30 

50 

100 

1500 

3000 

32.2 

* 

* 

407.0 

668.0 

of  object  speed.  Restricting  consideration  to  only  those  simulations 
resulting  in  sliding  transport,  Table  3  indicates  that  less  viscosity 
is  needed  at  low  overpressure  than  at  high,  since  n  =  300  g/s  appears  to 
be  adequate  for  a  cube  originally  at  30  psi;  at  Pgo  =  100  psi,  however, 
a  value  of  M  >1000  seems  required.  However,  another  approach  to  the  brak¬ 
ing  problem  appears  more  fruitful. 

Horizontal  Ground  Spring .  The  bouncing  mode  of  blast  transport  appears 
realistic  for  most  object-soil  interactions.  It  does  not  seem  at  all  likely, 
for  example,  that  the  grassy  surface  shown  in  the  photograph  contained  in 
Figure  1  would  allow  a  sliding  motion  of  even  a  heavy,  relatively  smooth 
object  such  as  the  concrete-filled  plywood  cube  observed  at  DIAL  PACK.  The 
torque,  therefore,  that  tends  to  tumble  a  blast-driven  object  appears  to  be 
a  desirable  feature  of  the  model.  To  be  able  to  keep  tumbling  transport,  yet 
get  enough  braking  on  ground  contact  to  simulate  experimental  results,  a 
horizontally-acting  spring  has  been  added  to  the  ground-object  interaction. 
Like  its  vertical  counterpart,  this  spring  is  double  valued  to  dissipate 
energy,  but  the  double  valuedness  also  serves  to  reduce  spinning. 

In  the  simplified  flowchart  of  DEBRIS  presented  in  Figure  10,  the 
box  marked  with  the  star  has  been  added  to  provide  notice  to  the  program 
when  one  or  more  corners  of  the  debris  fragment  first  come  into  ground 
contact,  at  which  time  logical  variable  GNDHIT  is  set  equal  to  .TRUE. 

At  the  same  time  GNDOFF  is  set  to  .FALSE,  so  that  on  the  next  time-cycle, 
GNDHIT  may  be  set  to  .FALSE.  These  logical  operations  permit  subroutine 
CORNER  to  record  the  coordinate  of  the  point  and  the  direction  of  hori¬ 
zontal  motion  of  most  recent  co  itact  and  to  measure  norizontal  spring 
compression  from  that  point.  Subroutine  CORNER  has  also  been  modified 
to  apply  the  horizontal  spring  forces  to  fragment  corners  below  ground. 

As  with  the  vertical  ground  spring,  horizontal  spring  constant  takes  a 
high  value  at  compression  and  a  low  value  at  release.  During  release, 
however,  should  the  corner  pass  through  the  point  of  most  recent  touch¬ 
down,  horizontal  spring  compression  may  begin  again,  the  direction  of 
horizontal  motion  noted,  and  the  high  value  of  spring  constant  be  appli¬ 
cable.  In  this  way,  considerable  energy  can  be  dissipated  upon  a  single 
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ground  contact — provided  the  fragment  remains  in  contact  and  does  not 
immediately  bounce  off.  The  direction  of  horizontal  spring  force  is 
always  pointed  toward  the  point  of  first  touchdown.  The  action  of  the 
ground  springs  is  sketched  schematically  in  Figure  11. 


Although  pure  friction  is  capable  of  inducing  bouncing  motion  in 
the  blast-transported  cube,  to  simulate  the  observed  transport  distance, 
the  coefficient  of  friction  must  vary  with  incident  overpressure  and 
increasing  the  coefficient  does  not  always  reduce  distance  travelled. 

From  the  limited  number  of  trials  performed  so  far,  it  appears  that 
using  horizontal  ground  springs  in  conjunction  with  the  vertical  the 
same  behavior  may  be  achieved  while  the  values  of  all  the  arbitrary 
parameters  are  held  constant.  These  parameters  must,  however,  be  changed 
with  the  soil.  They  probably  also  depend  on  the  nature  of  the  trans¬ 
ported  object.  In  the  work  to  date,  constant  values  of  friction  p  =  0.5 
and  of  viscosity  n  =  3000  g/s  have  been  left  in  the  calculations  through¬ 
out.  Spring  constants  have  been  varied,  however,  and  it  is  clear  that 
the  two  springs,  the  horizontal  and  the  vertical,  must  cooperate  in  pro¬ 
ducing  the  desired  motion.  Two  simulated  partial  trajectories  shown  in 

Figures  12  and  13  demonstrate  this.  Horizontal  constant,  k  ,  , ,  is  held 

ghl 

constant  in  the  two  calculations,  but  the  vertical  constant,  k  ,  ,  is 

gl 

changed  by  a  factor  of  ten.  The  stiffer  vertical  spring  causes  much 
higher  bounce  than  the  weaker  spring  and  the  greater  ground  contact  pres¬ 
ent  with  the  weaker  spring  causes  increased  loss  of  horizontal  speed. 

g 

When  k  ^  =0.35  x  10  dyne/cm, distance  moved  during  the  first  second  after 
blast  arrival  is  approximately  85  ft;  increasing  k  ^  by  a  factor  of  ten 
raises  the  transported  distance  in  the  same  time  interval  to  nearly  95  ft. 
The  loss  of  forward  momentum  might  be  caused  by  any  or  all  of  the  three 
mechanisms:  friction,  viscosity,  or  the  double-valued  horizontal  spring. 

Further  research  is  required  to  untangle  these  effects. 


With  the  parameter  values  unchanged  the  simulations  represented  in 
Figures  12  and  13  have  been  completed  and  two  aspects  of  these  calcula¬ 
tions  are  shown  in  Figures  14  and  15.  Incident  peak  overpressure  is  50 
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object  at  touchdown 


present  position 


touchdown 


Figure  11.  Schematic  Illustration  of  Ground  Reaction  Springs 
in  DEBRIS. 
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psi.  Figure  14  is  a  trajectory  of  the  center  of  mass  downrange  from 
which  it  is  seen  that  the  center  of  mass  never  rises  above  approximately 
1.5  in.  from  its  original  height.  The  small  squares  along  the  top  of 
the  figure  are  meant  to  represent  the  approximate  orientation  of  the 
cube  at  the  locations  shown;  the  simulated  cube  turns  completely  over 
once  before  reaching  the  top  of  its  first  bounce,  but  the  rate  of  tumbling 
declines  rapidly  after  two  bounces.  (This  is  in  contrast  to  behavior  in 
simulations  with  friction  or  viscosity  or  both  and  without  the  horizontal 
spring.)  Plotted  in  Figure  15  are  two  aspects  of  the  motion  illustrating 
braking  of  the  rotary  motion  and  of  the  translation.  The  peak  speeds  are 
reached  very  early  in  the  travel  under  the  influence  of  the  dynamic  pres¬ 
sure;  most  of  the  transport  is  during  the  braking  phase.  Maximum  rotary 
speed  is  approximately  9  rps  and  maximum  horizontal  speed  is  approximately 
38  ft/s.  These  values  depend  on  the  parameter  values  and  it  is  not  known 
from  the  experimental  data  whether  or  not  the  tumbling  motion  in  the 
simulation  is  an  accurate  representation  of  the  actual  motion.  The  sharp 
fluctuations  of  the  speeds  during  the  trajectory  are  not  real  and  the  dis¬ 
continuities  of  first  derivatives  are  due  mostly  to  the  lack  of  detailed 
computer  printout;  fluctuations  in  angular  and  translational  speeds  no 
doubt  actually  occur.  Also  no  attempt  has  been  made  as  yet  to  simulate 
the  fluctuation  in  transport  caused  by  variation  of  ground  surface  con¬ 
ditions.  In  the  DIAL  PACK  data,  for  example,  there  is  a  50  percent  vari¬ 
ation  in  observed  transport  at  50  psi. 

Figures  14  and  15  list  the  value  of  the  integration  tone  step  At. 

In  these  simulations,  when  the  transport  depends  critically  on  ground 
contact  during  bounces,  the  fineness  of  the  time  mesh  is  important.  The 
value  of  At  entered  in  the  figures  is  0.001  s,  but  during  the  time  any 
part  of  the  cube  is  in  contact  with  the  ground,  the  code  automatically  re¬ 
duces  the  preset  time  step  by  a  factor  of  ten  in  size  so  that  during 
most  of  the  motion  simulated  in  Figures  14  and  15  the  time  integration 
was  performed  with  steps  of  0.0001  s.  For  example,  increasing  the  nomi¬ 
nal  value  of  At  to  0.01  s  can  increase  calculated  transport  significantly 
in  some  over-the-ground  translations. 
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By  trial  and  error,  the  parameter  values 


k  =  0.20  x  10^  dyne/cm 
gl 

k  ^  =0.35  x  10^  dyne/cm 
p  =  0.5 

n  =  3000  g/s 


were  successful  in  simulating  distance  transported  in  the  range  of  in- 
cident  pressures  30  to  100  psi.  The  values  of  the  recovery  spring  con¬ 
stants  are  very  small  or  zero  and  do  not  seem  to  be  critical  to  the 
distance  travelled.  Distances  transported  in  simulations  with  these 
values  are  plotted  as  triangles  in  Figure  16  along  with  the  reported 
observations. 

At  least  in  this  pressure  range,  "tuning"  of  the  model  may  most 
easily  be  accomplished  by  adjusting  the  constant,  kgp  The  strong  influ¬ 
ence  of  the  vertical  compressive  spring,  k  ^ ,  on  braking  is  demonstrated 
by  the  distances  transported  reported  in  Table  4. 


Effect  of  Yield 

Because  the  DIAL  PACK  event  was  a  surface  explosion  of  approximately 
500  tons  of  high  explosive,  the  equivalent  nuclear  explosion  has  been 
simulated  by  airblast  arising  in  a  1-KT  surface  burst.  To  explore  the 
influence  of  yield  on  the  translation  model,  airblast  parameters  corre¬ 
sponding  to  a  1-MT  surface  burst  have  been  used  in  one  simulation  with 
the  result  that  total  travel  from  a  location  on  the  100-psi  contour  was 
2060  ft.  Simulation  parameter  values  were  the  optimum  values  noted  above. 


& 

At  15  psi,  the  diffraction  phase  is  relatively  more  effective  than  at 
higher  overpressures.  Allowing  for  diffraction  by  starting  the  cube 
at  t  =  1  ms  with  a  speed  of  1.84  ft/s  increases  displacement  by  approxi¬ 
mately  1  ft;  also  changing  kg2,  the  recovery  ground  spring  constant,  from 
0.35  x  1C)4  to  0.35  x  10^  dyne/cm  raises  transport  at  15  psi  to  6.26  ft, 
which  is  close  to  the  observed  average  7.5  ft.  See  appendix  for  a 

calculation  of  diffraction  momentum. 
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figure  16.  Simulated  and  Observed  Translations  of  Concrete  Cube  at  DIAL  PACK. 


This  value  of  transport  is  about  the  average  of  the  DIAL  PACK  observa¬ 
tions  scaled  from  1  KT  to  1  MT.  However,  to  save  computation  costs,  the 
time  siep  for  the  megaton  calculation  was  increased  to  0.01  s.  i’his 
change  probably  introduced  error. 

Comparison  with  Other  Data 

The  model  DEBRIS  can  simulate  the  translation  experiments  described 
by  Fletcher  and  Bowen" ">2 5  by  the  choice  of  the  initial  conditions  in  the 
calculation.  When  the  height  of  the  bottom  surface  of  the  dropped  object 
is  2.6  feet,  the  initial  speed  is  38  ft/s,  and  the  departure  time  later 
than  the  positive  dynamic  pressure  phase  duration,  the  calculated  motion 
with  the  optimum  simulation  parameters  noted  above  for  DIAL  PACK  does 
not  agree  with  the  described  behavior  of  the  stones  and  masonry  blocks 
dropped  from  the  moving  truck.  In  fact,  the  strong  horizontal  spring 
appears  to  throw  the  object  backward  toward  the  point  of  origin.  How¬ 
ever,  when  the  horizontal  spring  constant  is  reduced  by  two  orders  of 
magnitude  or  more,  the  computed  motion  is  more  realistic.  The  trajectory 
of  the  center  of  mass  according  to  a  simulation  of  the  drop  from  a  moving 
truck  for  the  concrete  cube  of  DIAL  PACK  is  drawn  in  Figure  17.  The  total 
distance  travelled  after  release  is  44. 6  ft.  According  to  Fletcher  and 
Bowen,22’23  the  distance  travelled  in  this  manner  is  independent  of  the 
n  >bs  of  the  object  dropped.  This  is  borne  out  by  DEBRIS,  when  a  cube 
identical  to  the  DIAL  PACK  cube,  but  having  a  density  of  0.08391  lb/in^ 
(which  is  the  density  of  normal  concrete),  is  simulated  with  DEBRIS.  In 
the  case  of  the  heavier  cube,  total  transport  is  44.0  ft.  The  trajectory 
(Figure  18)  is  again  realistic,  and,  in  fact,  quite  similar  to  that  ap¬ 
pearing  in  Figure  17  for  the  much  lighter  cube.  Whether  or  not  the  simi¬ 
larity  of  the  trajectories  is  accurate  is  not  known. 

In  any  case,  the  distance  and  time  of  transport  provided  bv  DEBRIS 
with  the  parameter  values  discovered  so  far  do  not  agree  with  the  pre¬ 
dictions  of  the  empirical  formulas  given  by  Fletcher  and  Bowen.22’23 
If  the  initial  speed,  v,  is  38  ft/s,  these  formulas  yield  distance,  s  = 
34.3  ft  and  time  to  stop,  t  =  1.907s.  In  the  case  of  the  heavy  cube 
DEBRIS  estimates  stopping  time  t  =  3.27  s  and  for  the  DIAL  PACK  density, 
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t  =  3.62  s.  Again  t  is  nearly  independent  of  mass,  as  forecast  by 
Fletcher  and  Bowen, 22,2^  but  the  quantity  differs  from  their  forecast. 

The  prediction  of  DEBRIS  as  to  distance  falls  within  one  standard  error 
of  estimate  of  the  value  predicted  by  Fletcher  and  Bowen's  formula;  the 
detailed  printout  shows  the  light  and  heavy  cubes  cease  downrange  progress 
at  2.7  s  and  3.3  s,  respectively. 
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PART  II 


BLAST  RESPONSE  AND  COLLAPSE  OF  BUILDINGS  (BRACOB) 

An  early  version  of  the  BRACOB  system  of  computer  programs  has  been 

2  5 

described  in  Rempel  and  Wiehle  and  a  then  current  listing  published  in 
Rempel.  The  system  accepts  information  about  the  floor  plan  in  terms 
of  a  user-defined  rectangular  coordinate  system,  engineering  parameters 
for  walls  and  airblast  specification  either  in  terms  of  a  classical  blast 
of  user-chosen  pressure  and  angle  of  impact,  or  in  terms  of  arbitrarily 
specified  pressure  histories.  The  program  itself  computes  loadings 
against  all  walls  as  functions  of  time,  response  of  walls  up  to  the 
point  of  collapse,  at  which  time  the  program  will  adjust  the  load¬ 
ing  on  other  walls  if  necessary.  Because  only  simplified  physical  prin¬ 
ciples  are  used,  i.e.,  structural  response  is  deduced  from  a  single  de¬ 
gree  of  freedom  model  and  airblast  loading  is  computed  with  the  help  of 
concepts  of  clearing  time  anu  room  filling,  the  computer  time  required 
for  a  complete  history  is  much  less  than  time  ordinarily  taken  for  finite 
element  or  hydrodynamic  calculations. 


The  single  degree  of  freedom  (SDOF)  model  provides  histories  of 
velocity  and  the  displacement  of  a  mass  equivalent  to  the  effective  mass 
of  the  wall,  which  is  assumed  to  bend  or  flex  under  an  airblast  load. 
Load-mass  factors  have  been  chosen  so  that  these  displacements  and 
speeds  correspond  to  observed  motions  of  central  points  in  the  responding 
walls. 2?j28  Further  evidence  of  the  applicability  of  the  SDOF  model 
to  estimates  of  deflection  are  expected  from  the  upcoming  MILL  RACE  event, 
particularly  in  regard  to  unreinforced  masonry.  Failure  of  an  unre¬ 
inforced  wall,  according  to  the  model,  takes  place  when  deflection  has 
reached  the  extent  that  the  wall  is  unable  to  support  overburden,  i.e., 


Particularly,  the  FEMA  Experiments  5401,  5402,  and  5403. 
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when  deflection  of  the  central  point  is  equal  to  the  thickness  of  the 
wall  itself.  At  this  point  BRACOB  ceases  to  compute  the  wall  movement. 

In  fact,  these  final  conditions  of  speed  and  deflection  become  initial 
conditions  for  the  code  DEBRIS;  provided  the  fracture  pattern  is  known 
as  well,  then  all  required  input  information  is  available  for  a  calcula¬ 
tion  of  the  continued  motion  of  the  wall  fragments  as  free  bodies. 

A  previous  study20  suggests  that  even  at  pressures  as  high  as  30 
psi,  unreinforced  masonry  fails  in  flexure,  i.e.,  that  shear  walls  pro¬ 
vide  support  of  some  degree  even  when  they  themselves  are  strongly  dis¬ 
turbed  by  the  blast.  (Again,  experiments  to  be  conducted  in  1981  at 
MILL  RACE  should  give  more  information  on  this  question.)  During  the 
present  research,  then,  exterior  unreinforced  walls  have  all  been  assumed 
to  fail  in  flexure,  demonstrating  typical  failure  patterns,  such  as  the 
classical  pattern  exhibited  in  Figure  19,  or  patterns  observed  during 
Shot  ENCORE.  0  ’20  The  ENCORE  walls  were  made  of  baked  clay  masonry 
units  within  steel  frames  and  were  exposed  in  the  range  from  4  to  9  psi. 
Portions  of  the  centers  of  these  walls  punched  out,  as  illustrated  in 
Figure  20. 

Published  reports  from  the  URS  blast  tunnel*  tests  of  masonry 
walls10’11'12  do  not  show  the  failure  patterns  of  Figures  19  and  20. 

In  these  experiments  the  test  walls  generally  filled  or  nearly  filled 
the  entire  cross  section  of  the  large  walk-in  gallery  that  was  used  as 
a  blast  tube;  furthermore,  support  conditions  at  wall  edges  were  often 
unusual.  Failure  in  walls  without  windows  appears  to  have  been  generally 
by  a  more  or  less  horizontal  crack  the  full  width  of  the  panel,  suggest¬ 
ing  that  vertical  edge  support  may  not  have  been  present,  or  not  effec¬ 
tive  . 


In  addition  to  wall  speed  and  time  of  failure,  the  DEBRIS  code  takes 
as  input  the  kind  of  hinging  before  collapse  or  assumes  none.  The  sub¬ 
sequent  trajectory  is  determined  in  large  measure  by  the  initial  condi¬ 
tions,  including  the  kind  of  hinging  or  its  absence. 


Located  at  Fort  Cronkhite,  Marin  County,  California. 


At  the  present  time  BRACOB  does  not  calculate  the  response  of  in¬ 
terior  walls,  but  it  will  treat  them  as  "collapsed"  as  of  a  certain 
preset  time  during  the  problem,  at  which  time  the  pressures  in  adjoining 
rooms  are  recalculated  by  means  of  thermodynamic  principles.  When  ex¬ 
terior  walls  collapse,  the  then  existing  outside  pressure  against  the 
wall  is  transferred  to  what  were  interior  walls  after  a  preset  time  (to 
allow  for  pressure  buildup  behind  the  collapsed  wall).  The  same  kind  of 
arbitrary  delay  is  imposed  on  the  room-filling  process  to  allow  for  (1) 
breakage  of  windows  and  doors,  and  (2)  interior  pressure  to  buildup 
throughout  the  room.  After  this  delay  on  room  filling,  BRACOB  calculates 
an  average  interior  pressure  resulting  from  the  air  streaming  in  from  out¬ 
side  high  pressure  regions. 

Because  reinforced  wall  panels  seldom  support  overburden  (except  in 
earthquake  zones),  the  "failure"  criterion  for  these  walls  is  different 
from  that  for  unreinforced  load-bearing  walls.  Failure  of  reinforced 
concrete  walls,  for  example,  is  defined  to  occur  when  a  certain  degree 
of  plastic  deformation  of  the  reinforcement  has  occurred.  Exactly  at 
what  point  in  their  response  to  air  blast  such  panels  become  "debris"  is 
not  known  at  this  time,  but  FEMA  Experiment  5402  at  MILL  RACE  should 
provide  some  insight  on  this  question. 

Because  of  the  greater  degree  of  understanding  of  unreinforced  wall 
behavior  after  collapse,  attention  during  this  research  has  focussed  on 
unreinforced  masonry  walls. 

Improvements  in  BRACOB 

In  preparation  for  systematic  studies  of  debris  formation  and  scatter 
in  cities,  BRACOB  has  been  enhanced  to  handle  all  angles  of  impact,  0  to 
2tt;  necessary  restructuring  of  the  data  base  when  incidence  is  from  the 
rear  is  now  accomplished  automatically. 

Better  use  of  the  powerful  computing  system  at  the  Stanford  Center 
for  Information  Processing  is  now  made  possible  by  a  system  of  independent 
compilation  of  subroutines.  Data  bases  also — for  specifying  airblast  and 
for  describing  the  building — have  been  individualized,  as  have  output 
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for  various  kinds  of  output  information.  Although  computation  remains 
in  the  batch  mode  to  take  advantage  of  reduced  costs,  operation  of  the 
system,  including  disposition  of  files,  may  be  interactive. 

Multiply-connected  spaces  in  floor  plans  formerly  caused  an  error 
termination  from  BRACOB;  this  shortcoming  has  been  corrected.  Storage 
specifications  have  been  changed  to  permit  calculation  of  blast  response 
of  larger  buildings  than  before. 

Another  enhancement  to  facilitate  application  of  BRACOB  to  large 
urban  buildings  is  the  creation  of  "loading  multipliers"  that  enable  the 
program  to  calculate  clearing  times  from  geometric  data  pertaining  to  a 
single  floor  of  a  multiple-story  building.  Clearing  time,  tc,  for  re¬ 
flected  pressure  against  a  facade  of  any  building  may  be  found  from  the 

O  Q 

formula  • 


n  ( A  . .  -  A  ) 
wall  opng 

c  aP 

where  A  , ,  is  the  total  area  of  the  facade  and  A  is  the  area  of 
wal 1  opng 

openings  in  the  wall,  a  =  sound  speed  in  the  reflected  zone  in  front  of 
the  facade,  and  P  =  the  total  perimeter  in  the  facade  from  which  clearing 
relief  waves  may  originate.  (Window  and  facade  edges  are  treated  as 
equal  despite  the  fact  that  room  interiors  may  contain  a  higher  pressure 
than  found  at  the  edge  of  the  building.) 

The  quantity  n  is  3  if  incident  overpressure  is  below  12  psi,  and 
n  =  4  if  incident  overpressure  is  >12  and  <25  psi.  There  are  three  mul¬ 
tipliers:  one  for  wall  area,  another  for  window  or  opening  area,  and 

a  third  for  perimeter.  Were  all  multipliers  unity,  the  calculation  of 
clearing  time  would  be  correct  for  a  single-story  building.  A  six-story 
building  with  six  identical  floors  would  require  all  multipliers  to  be 
equal  to  6  because  the  roof  is  claimed  as  free  perimeter  only  once,  i.e., 
the  perimeter  multiplier  is  applied  only  to  the  length  of  vertical  edges. 
In  the  case  of  oblique  incidence,  the  vertical  facade  edge  is  a  free 
perimeter  for  only  one  of  the  two  adjoining  facades. 
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Certainly  one  of  the  major  enhancements  accomplished  during  this 
period  of  research  has  been  the  incorporation  into  the  subroutine  RESIST 
(which  determines  the  parameters  of  the  resistance  function  for  each  ex¬ 
terior  wall)  of  provision  for  Support  Conditions  9  and  10. 28  Walls  so 
supported  are  unreinforced  and  arch  in  a  flexible  or  rigid  frame,  a  con¬ 
dition  common  among  large  urban  buildings  in  the  United  States  outside 
California  or  other  earthquake-prone  areas.  Support  Condition  9  refers 
to  one-way  arching  in  a  frame  and  Condition  10  describes  two-way  arching. 
The  arching  code  was  originally  written  to  calculate  the  blast  response 
of  half-timbered  walls,  i.e.,  bricks  inside  a  wooden  frame,  in  which 
case  the  yielding  of  the  frame  was  important  to  the  result. 

In  steel  frame  or  even  reinforced  concrete  frame  buildings,  yielding  of 
the  frame  is  not  usually  significant.  For  the  present  research  the 
frame  has  been  considered  rigid. 

LANDIS  HOSPITAL 

The  Henry  R.  Landis  State  Hospital,  Philadelphia,  Pennsylvania,  is 
one  of  the  buildings  in  the  National  Shelter  Survey  inventory.  An  ele¬ 
vation  view  of  one  wing,  taken  from  blueprints  received  from  Research 
Triangle  Institute  appears  in  Figure  21  and  a  simplified  typical  floor 
plan  in  Figure  22  showing  the  orthogonal  grid  required  by  BRACOB  as 
the  basis  for  geometric  input  data.  The  computer  program  BRACOB  has  pro¬ 
duced  the  schematic  drawing  in  Figure  23.  The  drawing  corresponds  to  a 
time  0.01  s  after  blast  arrival;  the  current  location  of  the  blast 
front  is  shown  by  the  dashed  line  across  the  first  tier  of  rooms.  Open¬ 
ings  are  represented  by  short  lines  that  parallel  the  walls.  As  can  be 
seen  from  the  figure,  EPJ\C0B  now  accepts  cases  in  which  rooms  are  en¬ 
tirely  surrounded  by  other  rooms;  that  is,  artificial  walls  are  not 
necessary  to  connect  the  inner  core  of  rooms  to  the  surrounding  rooms. 

In  Figure  23>  the  blast  front  parallels  the  long  axis  of  the  floor  and 
travels  from  the  bottom  to  the  top  of  the  page. 

As  can  be  seen  from  the  structural  detail  of  the  outer  walls  repro¬ 
duced  in  Figure  24,  there  are  two  tiers  of  masonry,  each  arching  in  a 
nearly  rigid  frame.  The  outer  tier  of  red  brick  is  framed  top  and  bottom 
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Landis  State  Hospital,  Philadelphia,  South  and  West  Elevation 
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Figure  22.  Orthogonal  Grid  for  Landis  Hospital. 
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Figure  24.  Outside  Wall  Detail,  Landis  Hospital 


bv  metal  angles  embedded  in  the  reinforced  concrete  spandrel  beam;  the 
inner  tier  of  concrete  masonry  units  arches  at  top  and  bottom  within  the 
R/C  frame  itself.  There  are  vertical  frame  members  only  at  every  other 
cross  partition  along  the  "front"  outside  wall,  i.e.,  the  wall  at  the 
bottom  of  the  page  in  Figures  22  and  23.  The  nonstructural  interior 
partitions  offer  strong  initial  support  to  the  outer  wall,  but  arching 
is  not  assumed  there;  hence,  the  exterior  front  walls  are  treated  in  the 
analysis  of  blast  response  as  "one-way  arching"  in  a  rigid  frame.  Figure 
24  shows  metal  ties  24  in.  on  center  vertically  between  the  two  tiers  of 
the  outer  walls;  this  is  not  regarded  as  enough  to  transfer  shear  forces, 
and  the  two  wythes  are  treated  analytically  as  independent  or  additive; 
that  is,  the  resistance  to  deflection  is  the  sum  of  the  two  acting  in¬ 
dependently.  Because  both  wythes  are  the  same  thickness,  "failure"  is 
equated  to  a  deflection  of  the  center  of  the  wall  equal  to  a  brick  width, 
i.e. ,  3.625  in. 

The  crushing  or  compressive  strength  for  the  outside  walls,  which, 
when  arching  within  a  rigid  frame,  determines  the  ultimate  resistance  of 
the  wall  to  flexure,  has  been  chosen  to  be  that  of  type  S  field-prepared 
mortar,  i.e.,  1600  psi.28  (This  value  is  also  approximately  that  for 
the  concrete  masonry  unit  itself,  but  the  clay  brick  is  stronger  in  com¬ 
pression  than  mortar.)  Except  for  the  variations  in  window  sizes,  all 
the  outside  walls  are  the  same  structurally.  Because  most  of  the  windows 
occupy  such  a  relatively  large  portion  of  the  wall  areas,  the  walls 
themselves  have  been  treated  as  one-way,  i.e.,  the  support  given  to  the 
walls  by  adjacent  edges  is  not  considered.  The  resistance  function  of 
an  outside  wall  of  the  Landis  building  it  computed  as  the  sum  of  the 
individual  resistances  of  two  identical  tiers,  each  3-5/8  in.  thick, 

1600  psi  in  compressive  strength,  modulus  of  elasticity  1  x  10^  psi, 
and  density  equal  to  the  sum  of  brick  and  concrete  masonry  units,  i.e., 
109.2  pcf,  plus  48.3  pcf  or  157.5  pcf.  (Taking  wall  thickness  equal  to 
"ha*:  of  a  single  tier  ensures  that  the  collapse  criterion  will  be  correct 
during  calculation  and  adding  the  densities  of  the  two  tiers  then  makes 
the  inertial  mass  correct.  Were  the  two  tiers  capable  of  bending  as  a 
unit,  i.e.,  were  the  metal  ties  capable  of  transmitting  shear  forces 


between  tiers,  then  the  total  resistance  of  the  wall  to  flexure  would 
be  considerably  higher  than  merely  the  sum  of  the  resistance  of  each  wal 
singly.  Using  the  sum  is  equal  to  considering  the  two  tiers  as  SDOF  sys 
terns  in  parallel,  i.e.,  that  the  ties  maintain  the  distance  between  the 
tiers  fixed,  so  that  the  deflection  of  the  two  SDOF  systems  are  equal.) 

Inside  Walls 

To  some  extent  the  behavior  of  the  outside  walls  depends  on  the  re¬ 
sponse  of  the  interior  walls.  Should  the  inside  walls  remain  in  place 
long  enough  to  hold  the  backfilling  pressure  against  the  inner  surfaces 
of  the  outside  walls  for  an  appreciable  time,  then  the  outside  walls 
will  appear  stronger  than  if  the  interior  walls  disappear  immediately 
and  open  up  the  whole  interior  space  to  blast  pressure.  In  its  present 
form,  BRACOB  will  remove  the  interior  walls  at  any  preset  time  after 
blast  arrival  at  the  wall,  but  it  cannot  now  analyze  the  inside  walls 
for  response  to  airblast.  In  the  case  of  the  Landis  Hospital,  it  does 
not  appear  that  the  interior  partitions  have  significant  structural 
strength,  except  as  shear  walls  (i.e.,  in  resisting  loads  applied  in  the 
direction  of  their  long  axes)  and  except  insofar  as  they  have  mass  to 
be  moved  out  of  the  path  of  the  airblast.  Therefore,  building  response 
has  been  analyzed  two  different  ways:  first,  all  interior  walls  lemain 
intact  throughout  the  calculation,  and  second,  interior  walls  are  re¬ 
moved  Instantly  at  blast  arrival.  The  corresponding  difference  in  re¬ 
sponse  of  tiic  outside  walls  is  measureable,  but  not  large.  Although  no 
attempt  has  been  made  to  pinpoint  it  exactly,  the  "incipient  collapse 
overpressure"  under  the  first  assumption  is  higher  than  under  the  second 
When  interior  walls  are  kept  in  place,  the  front  walls  reach  failure  de¬ 
flection  at  a  loading  applied  by  a  1-psi  head-on  blast  between  48  and 
50  ms  after  blast  arrival;  furthermore,  they  are  at  that  time  being 
pushed  weakly  outward  by  the  slightly  higher  inside  pressure;  on  the 
other  hand,  without  the  presence  of  the  interior  walls,  but  with  the 
same  incident  blast  conditions,  the  front  walls  fail  at  39  ms,  when  they 
are  still  accelerating  inwaidlv.  Departure  speeds  of  the  center  of  mass 
of  the  front  wall  fragments  vary  somewhat  from  room  to  room,  depending 
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on  wall  span,  room  volume,  and  window  area,  but  under  the  first  assump¬ 
tion  these  speeds  vary  from  6.04  to  7.15  fps  (1.84  to  2.18  m/s).  Under 
the  second  assumption  about  the  behavior  of  the  interior  walls,  the 
variation  is  from  11.45  to  12.4  fps  (3.49  to  3.79  m/s). 

Because  each  window  in  the  front  wall  occupies  such  a  large  portion 
of  its  wall,  the  failure  pattern  will  be  close  to  that  sketched  in  Fig¬ 
ure  25.  The  four  fragments  will  most  likely  be  temporarily  "hinged"  at 

their  outer  edges,  that  is,  will  tend  to  rotate  initially  about  the  edge. 

Under  such  circumstances,  their  trajectories  will  take  them  into  collision 
with  a  ceiling,  sidewall,  or  floor,  all  of  which  structural  elements  are 
very  likely  to  be  still  in  place  at  the  time  of  collision  with  the  front 

wall  debris.  The  effect  of  this  debris  mass  on  these  surfaces  is  not 

clear  at  this  time.  The  impacts  represent  transient  loadings  in  excess 
of  design  capacities. 

Collapse  Scenarios 

The  collapse  of  walls  and  creation  of  wall  debris  in  the  Landis 
Hospital  are  described  by  means  of  the  foregoing  techniques. 

The  Landis  Hospital  is  regarded  as  a  typical  large  urban  building 
found  outside  the  earthquake  zones  of  the  United  States.  With  small 
changes  in  details,  the  scenarios  presented  here  are  expected  to  be 
similar  to  the  airblast  responses  of  nearly  all  buildings  of  this  class. 

Two  observations  about  the  Landis  Hospital  are  appropriate,  however. 
First,  the  exterior  walls  are  relatively  weak.  Even  the  assumption  of 
full  arching  in  a  rigid  frame  does  not  raise  their  strength  significantly 
because  each  tier  is  so  thin.  A  4-in.  thick,  unreinforced  masonry  wall 
does  not  resist  airblast;  even  adding  two  of  them  together  does  not  cre¬ 
ate  a  good  blast  shelter.  Second,  the  relatively  large  windows  in  the 
south  and  north  elevations*  permit  a  rapid  buildup  of  interior  pressure 
in  the  first  tier  of  rooms,  but  the  relatively  large  open  spaces  on  the 


Airblast  approaching  from  either  the  south  or  north  will  have  almost 
the  same  effect  on  Landis  Hospital. 
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far  side  of  the  first  tier  of  rooms  means  that  net  pressure  on  some  of 
the  interior  walls  becomes  large  enough  to  break  their  connections  to 
the  frame  and  convert  them  into  debris,  although  the  interior  walls  them¬ 
selves  may  remain  intact  for  some  time  (as  described  in  Gabrielsen  et  al.11) 
The  interior  walls  are  similar  both  in  method  of  anchorage  and  in  construc¬ 
tion  to  interior  walls  exposed  to  low-pressure  airblast  in  the  'JRS  shock 
tunnel  during  the  early  1970s.11  Figure  26  is  taken  from  the  construction 
drawings  for  Landis  Hospital  to  illustrate  the  interior  walls. 

The  most  nearly  comparable  panels  tested  in  the  URS  tunnel  were  made 
of  metal  studs  covered  with  4-in.  sheetrock,  fastened  to  the  rock  tunnel 
wall  with  explosively  driven  studs.  These  examples  were  tested  at  only 
two  peak  incident  overpressures,  1.7  and  3.5  psi.  The  results  were  re¬ 
ported  in  the  form  of  photographs  before  and  after  and  measured  displace¬ 
ment  histories.  These  walls  apparently  travelled  outward  from  their 
original  positions  as  single  units,  at  least  for  several  feet.  Displace¬ 
ment  histories  for  the  two  pressures  are  shown  in  Figure  27  (Gabrielsen 
et  al.11).  Figure  27  also  contains  two  conjectural  histories  for  1  and 
30  psi,  calculated  from  the  reported  data  by  simple  interpolation  (a  pro¬ 
cedure  that  essentially  neglects  the  strength  of  the  fasteners).  The 
conjectural  curves  in  Figure  27  will  be  used  later  in  the  analysis. 

Certain  major  events  and  conditions  are  predicted  as  a  nuclear- 
derived  airblast  impinges  directly  against  the  south  elevation,  reflects 
and  diffracts  around  the  structure,  and  finally  sweeps  through  the  in¬ 
terior.  For  the  purposes  of  this  discussion,  the  airblast  originates  in 
the  explosion  of  1  MT  on  the  surface.  Incident  pressures  are  restricted 
to  1  and  30  psi,  providing  two  extremes  in  response. 

The  influence  of  backloading  on  the  collapse  of  the  outside  walls 
of  the  south  elevation  is  clarified  by  Figures  28  and  29  .  In  both  fig¬ 
ures  the  exterior-reflected  pressure  as  a  function  of  time  is  shown  by 
Curve  A.  In  Figure  28  the  interior  pressure  within  the  first  tier  of 
rooms  (e.g.,  Room  Number  2  in  Figure  23),  shown  by  Curve  B,  rises  to  its 
peak  in  approximately  30  ms  at  which  time  inside  room  pressure  sur¬ 
passes  outside  pressure  slightly.  However,  enough  momentum  is  imparted 
to  the  wall  by  the  short-lasting  blast  impact  that  the  wall  fails  later 
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Figure  26.  Interior  Wall  Detail,  Landis  Hospital. 


Interior  and  Exterior  Pressures,  Landis  Hospital,  All  Walls  Present 


wall  it 4  fails  Inward  (front  tier  rooms  absent) 
v  =  7.81  fps 


Figure  29.  Interior  and  Exterior  Pressures,  Landis  Hospital  Front  Rooms  Removed, 


at  approximately  72  ms  with  an  inward  speed  of  4.52  fps.  In  Figure  29 
the  interior  pressure  has  been  calculated  for  the  case  of  no  interior 
walls  in  the  front  tier  of  rooms.  The  initial  floor  plan  for  this  cal¬ 
culation  is  represented  in  Figure  30.  Under  these  conditions  the  back- 
loading  is  still  rising  as  late  as  55  ms,  and  the  high  net  pressure  (i.e., 
Curve  A  less  Curve  B)  collapses  the  outside  walls  in  the  south  facade 
at  approximately  52  ms  and  imparts  speeds  of  7.81  fps.  Probably  neither 
of  the  two  Curves  B  is  correct  because  of  the  presence  of  the  interior 
walls  in  the  front  rooms  for  a  portion  of  the  time  during  exterior  load¬ 
ing.  It  is  conjectured,  based  partly  on  the  work  reported  in  Cabrielsen 
et  al.,11  that  the  interior  partitions  break  their  anchorages  at  least 
as  soon  as  the  pressure  differential  across  them  reaches  1  psi  and  that 
at  that  moment  equalization  of  air  pressure  between  the  front  tier  of 
rooms  and  the  open  spaces  behind  them  begins.  In  Figure  28  ,  the  com¬ 
puted  pressure  in  this  open  space  (designated  as  Room  21  in  Figure  23) 
calculated  by  assuming  the  presence  of  all  interior  walls  is  shown  by 
Curve  C.  Clearly,  backfilling  into  Room  21  does  not  occur  fast  enough 
to  prevent  failure  of  most  of  the  interior  walls  in  the  front  tier  of 
rooms  (Numbers  1  through  20  in  Figure  23). 

Curve  D  in  Figure  28  represents  the  backloading  against  the  front 
exterior  walls  under  the  assumptions  that  at  approximately  20  ms  the  ex¬ 
posed  interior  walls  start  to  move  out  of  the  way  of  the  air! last  enter¬ 
ing  through  the  front  windows,  and  that  at  30  ms  the  pressure  is  equalized 
throughout  the  area  indicated  as  Room  1  in  Figure  30.  (The  shear  walls 
supporting  the  13  individual  front  walls  in  the  south  elevation  are  re¬ 
placed  by  pilasters  in  the  calculations  based  on  Figure  30.)  Using  a  net 
load  equal  to  Curve  A  less  Curve  D,  response  of  the  front  walls  can  be 
calculated.  It  is  clear  the  calculation  will  predict  the  conditions 
of  exterior  front  wall  collapse  between  those  given  in  Figures  28  and  29. 
In  estimating  debris  formation  and  distribution,  an  average  will  be  used, 
i.e.,  failure  at  62  ms  and  fragment  speed*  equal  to  6.16  fps. 


The  speed  of 


the  center  of  mass  of 


the  fragments. 
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Figure  30-  Special  Initial  Floor  Plan,  Landis  Hospital. 


The  rooms  in  the  central  core  (Numbers  2  through  9  in  Figure  30) 
will  most  likely  disappear  soon  after  the  blast  reaches  them,  unless 
there  are  large  openings  in  the  walls,  in  which  case  they  may  survive 
the  diffraction  phase  to  perish  in  the  later  drag  phase. 

No  loading  of  an  interior  wall  can  begin  before  the  blast  wave  can 
reach  it;  hence,  the  earliest  any  load  can  be  expected  against  the  in¬ 
terior  walls  defining  the  rear  tier  of  rooms  (along  Line  I  =  11  in  Fig¬ 
ure  22)  is  45  ms.  However,  after  this  time,  the  pressure  builds  up  much 
faster  on  the  south  side  of  this  line  of  partitions  (facing  the  oncoming 
blast)  than  it  does  on  the  north  side,  facing  downstream.  This  is  shown 
in  Figure  29  by  the  two  Curves  E  and  F,  which  represent  the  outside  and 
inside  pressures,  respectively,  in  the  rear  tier  of  rooms  (Numbers  10 
through  21  in  Figure  30  ) .  Evidently  the  interior  and  exterior  pressures 
remain  nearly  in  equilibrium  because  of  the  relatively  slow  outside 
buildup  and  the  large  windows  in  the  north  facade.  It  is  presumed  that 
the  interior  partitions  separating  the  rear  tier  of  rooms  from  the  rest 
of  the  floor  will  fail  at  a  time  between  60  and  70  ms,  and  pressure 
within  those  rooms  will  rise  to  the  values  given  in  Figure  29  by  Curve  G. 
Under  these  conditions,  the  outside  walls  of  the  north  facade  are  pre¬ 
dicted  to  collapse  outward  at  time  equal  to  170  ms  with  a  departure  speed 
of  approximately  3.4  fps.  The  floor  is  thus  completely  opened  to  air- 
blast  at  a  time  when  the  free-field  static  pressure  still  has  95  percent 
of  its  peak  value.  Further  downrange  translation  of  wall  debris  can  be 
expected . 

Because  of  the  temporary  shielding  of  the  interior  offered  by  the 
front  tier  of  rooms,  under  certain  conditions  some  side  walls  may  fail 
as  well  as  the  forward  and  backward  facing  walls.  For  example,  Wall 
Number  5  (Figure  30) ,  will  not  be  backloaded  until  the  partition  line  I  = 
11  ceases  to  block  blast  pressure.  Calculations  show  that  Wall  Number  5, 
without  backloading,  fails  inward  at  103  ms.  However,  as  noted  above. 
Wall  Line  11  may  well  become  ineffective  before  103  ms,  in  which  case 
the  motion  of  Outside  Wall  Number  5  will  be  slowed,  if  not  reversed.  At 
1  psi  incident  pressure,  all  exterior  sidewalls  either  do  not  fail  during 

*  Assuming  inside  pressure  rises  due  to  filling  through  windows  in  the 
north  facade. 
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the  period  covered  by  calculations  to  date,  or  do  so  marginally,  as 
Wall  Number  5. 

The  comparable  analysis  for  an  incident  overpressure  of  30  psi 
presents  a  much  faster  mo\ ing  version  of  the  same  picture  as  the  analysis 
for  PSQ=  1  psi.  The  pressures  responsible  for  the  destruction  are 
depicted  in  Figure  31.  Curves  A  and  F  represent  the  two  exterior  pres¬ 
sure  histories,  A  against  the  south  facade  and  F  against  the  north 
facade.  The  remainder  of  the  curves  in  Figure  31  show  interior 
pressures  calculated  under  various  assumptions.  Curve  B  is  the  backload¬ 
ing  against  the  south  facade,  assuming  all  the  interior  walls  are  in 
place  and  act  as  effective  barriers  to  rapid  pressure  rise.  The  insig¬ 
nificant  level  of  the  backloading  pressure  contrasted  to  the  exterior 
pressure  means  almost  instantaneous  collapse  of  the  south  facade;  in  the 
south  facade,  calculation  shows  failure  6  ms  after  blast  arrival  and  a 
speed  of  departure  of  the  center  of  mass  of  the  wall  fragments  equal 
to  111  fps.  The  interior  walls  defining  the  front  tier  of  rooms  dis¬ 
appear  within  a  few  milliseconds  because  they  are  significantly  loaded 
by  simple  room  filling  through  the  windows  in  the  south  facade;  the  ad¬ 
vent  of  the  high  reflected  pressure  against  the  south  facade  upon  failure 
of  that  facade  at  approximately  16  ms  will  shear  the  connections  of  the 
interior  panels  within  much  less  time  than  the  same  task  required  during 
the  URS  tunnel  tests.11  The  conjectural  displacement  history  correspond¬ 
ing  to  p  =  30  psi  in  Figure  27  suggests  a  separation  from  the  frame 
equal  to  1  ft  within  10  ms  of  blast  arrival  at  the  partition.  Because 
the  exterior  blast  front  reaches  the  grid  line  I  =  11  at  approximately 
25  ms,  the  surge  of  high  pressure  resulting  from  the  failure  of  the 
south  facade  may  be  expected  at  the  interior  partition  line  I  =  11  about 
10  ms  later  or  at  time  equal  to  35  ms.  This  is  before  any  backfilling 
of  the  north  tier  of  rooms  can  begin  (shown  by  Curve  D  in  Figure  31)  and 
after  average  interior  pressure  in  the  remainder  of  the  floor  has  reached 
approximately  10  psi  (Curve  E) .  Again,  the  interior  partition  line  at 
I  =  11  will  cease  acting  as  a  pressure  barrier  within  10  ms,  admitting 
high  pressure  at  that  time  to  the  rear  tier  of  rooms.  This  high  pres¬ 
sure  is  between  values  given  by  Curves  A  and  C,  that  is,  an  order  of 
magnitude  higher  than  the  backfilling  pressure  entering  through  the 

windows  in  the  north  facade.  The  exterior  walls  of  the  north  facade  then 
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begin  their  response  to  the  30-psi  blast  by  moving  inward  as  the  blast 
engulfs  the  building,  but  almost  immediately  their  acceleration  is 
reversed  because  of  the  rapid  buildup  of  interior  pressure  and  these 
exterior  walls  ultimately  collapse  outward.  Current  estimates  place 
the  time  of  failure  in  the  north  facade  at  62  ms  and  outward  speed  of 
the  fragments  at  8  fps. 

North  and  south  exterior  wall  failure  data  are  summarized  in  Table  3. 

The  behavior  of  the  shear  walls  at  30  psi  is  not  so  clear.  Under 
the  influences  of  ordinary  room  filling  through  the  windows  in  the  south 
facade  and  blast  loading  as  the  front  passes  down  its  length,  the  two 
first-loaded  shear  walls  (Numbers  3  and  117  in  Figure  23  )  move  inward 
and  fail  at  23  and  17  ms,  respectively,  and  with  speeds  of  27  and  40  fps, 
respectively.  However,  the  front  wall  opens  up  after  6  ms  and  interior 
pressure  in  the  front  tier  of  rooms  rises  more  rapidly  than  shown  in 
Curve  B;  in  fact,  the  effective  pressure  increase  is  probably  similar  to 
that  suggested  by  the  two  Curves  G  in  Figure  31.  However,  even  the 
steeper  of  the  two  Curves  G  will  not  reverse  the  ultimate  inward  failure 
of  the  shear  walls  in  the  fron  tier  of  rooms,  although  the  wall  speed  at 
failure  will  depend  on  the  exact  assumption  regarding  interior  pressure. 

Calculation  suggests  that  the  behavior  of  the  front  outside  shear 

walls  is  typical  of  their  class  at  p  =30  psi.  Unless  it  becomes  evi- 

so 

dent  from  further  experimentation,  interior  blast  loading  does  not  appear 
rapid  enough  to  halt  the  strong  inward  motion  imparted  to  outside  shear 
walls  by  the  initial  blast  impact.  Walls  Number  1  and  109  fail  inwardly, 
for  example,  at  32  ms  with  computed  speeds  near  45  fps.  Such  a  compu¬ 
tation  admittedly  is  based  on  ordinary  room-filling  and  does  not  allow 
for  blast  removal  of  exterior  walls;  but  the  time  of  failure  is  so  early 
and  the  speed  so  high  that  it  is  not  likely  to  be  reversed. 

Inward  failure  of  shear  walls  does  not  necessarily  imply  inside 

20 

deposition  of  shear  wall  debris,  a  subject  that  will  be  discussed  later. 
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Table  5 


FAILURE  OF  EXTERIOR  WALLS,  LANDIS  HOSPITAL 


Yield  =  1MT 

HOB  =  0 

Peak  Free-field  ^ 

\  South 

North 

Pressure  (psi) 

\  Elevation 

Elevation 

Inward 

Outward 

1 

v^  =  6.16  fps 

3. A  fps 

t  =  62  ms 

c 

170.0  ms 

Inward 

Outward 

v^  *  111.  fps 

8.0  fps 

30 

t  =  6  ms 

c 

62 .  ms 

Vj  =  speed  of  departure  of  fragments 

t  =  time  of  structural  failure 
c 
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Fragment  Distribution,  pgo  =  1  psi 

After  the  failure  of  the  exterior  walls  in  the  south  elevation  there 
is  continued  flow  into  the  building  as  the  blast  finds  its  way  deeper  in¬ 
to  the  interior.  Such  flow  helps  carry  the  fragments  of  the  failed  walls 
but  the  major  contributor  to  their  momentum  at  this  time  is  still  the  pre 
vious  impact  of  the  air  blast. 

The  flow  stream  has  been  established  around  the  structure.  There 
are  only  weak  or  short-lasting  pressure  differentials  within  the  building 
itself.  However  at  some  time  a  flow  path  is  opened  up  through  the  build 
ing;  this  will  occur  after  the  failure  of  the  north  exterior  walls  at  170 
ms  makes  a  connection  between  the  high  pressure  on  the  south  and  the  re¬ 
latively  low  pressure  on  the  north.  Gradually  a  flow  pattern  through  the 
building  will  be  established.  This  extremely  complicated  process  will  be 
represented  in  these  simplified  calculations  by  a  hiatus  of  zero  flow 
immediately  after  the  failure  of  the  southern  external  facade  followed 
by  free-field  flow  beginning  after  the  failure  of  the  northern  facade 
is  communicated  to  the  upstream  facade  of  the  building  by  a  sound  signal. 
At  1  psi  under  standard  conditions  sound  speed  is  1150  fps;  the  building 
is  76  feet  wide;  so  that  the  delay  in  establishing  flow  after  failure  of 
the  northern  facade  will  be  assumed  to  be  66  ms.  Thus  the  front  wall 
fragments  will  travel  for  170  +  66  -  62  =  174  ms  without  acquiring  momen¬ 
tum  from  air  flow.  Fragments  at  the  northern  facade  will  do  the  same  for 
twice  66  ms  or  132  ms  after  failure.  Upon  expiration  of  the  hiatus 
classical  air  blast  flow  appropriate  to  the  adjacent  free-field  environ¬ 
ment  will  be  assumed  to  translate  the  fragments. 

Calculation  with  DEBRIS  shows  that  during  the  hiatus  the  south  wall 
fragments  translate  approximately  0.6  ft  downstream,  fall  approximately 
0.5  ft  toward  the  floor  and  rotate  40  degrees.  Since  the  fragments  will 
soon  strike  interior  structures  such  as  the  floor,  shear  walls  and  ceil¬ 
ing,  the  precise  orientation  at  any  time  is  not  important,  and  the  trans¬ 
lation  distances  during  the  hiatus  are  clearly  negligibli  .  In  calcul¬ 
ation  of  the  further  displacement  of  the  fragments  from  the  front  wall 
the  greatest  uncertainty  lies  in  the  choice  of  parameters  characterizing 
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the  concrete  floor  for  the  purposes  of  the  DEBRIS  code.  No  experimental 

data  are  known  from  which  to  deduce  values  for  the  spring  constants, 

viscosity  and  friction.  Since  the  horizontal  ground  spring  seems  to  be 

most  useful  in  representing  rough  terrain  rather  than  smooth  surfaces, 

it  will  not  be  incorporated  into  the  calculations  of  front  wall  debris 

transport.  The  vertical  spring  should  presumably  be  stiff;  thus  k  will 

9 

be  set  equal  to  0.35  x  10  dyne/cm.  Friction  and  viscosity  will  be  un¬ 
changed  from  the  values  y  =  0.5  and  p  =  3000  g/s.  Initially,  no  con¬ 
sideration  will  be  taken  of  the  presence  of  other  debris  within  the 
building  nor  of  the  possible  further  breakup  of  the  four  individual  frag¬ 
ments  into  which  the  wall  is  assumed  to  have  broken  by  the  blast  impact. 
(The  fragmentation  pattern  is  that  of  Figure  25.) 

Detailed  results  of  calculations  with  DEBRIS  show  that,  after  re¬ 
sumption  of  classical  air  blast  flow  through  the  building,  the  south  ex¬ 
terior  wall  fragments  fall  into  the  concrete  floor  A  to  7  ft  inside  the 
wall  depending  on  origin*  and  that  none  of  these  fragments  under  the  pre¬ 
sent  assumptions  advances  further  than  2  ft  from  the  point  of  impact. 

The  fragments  of  the  northern  exterior  wall  fall  directly  into  the  ground 
outside  the  building  and  even  under  the  assumptions  used  for  the  simu¬ 
lation  of  the  concrete  surface  these  fragments  do  not  bounce  more  than  a 
few  inches.  This  behavior  appears  to  be  due  to  the  simultaneous  absence 
of  strong  rotation  and  strong  dynamic  pressure  at  the  time  of  impact. 

(A  fragment  rotating  into  a  surface  has,  as  noted  before,  the  opportunity 
to  be  propelled  upward  by  the  compressive  ground  spring  after  which  it 
may  be  caught  and  rotated  further  by  blast  wind,  if  any  exists  at  that 
time.)  Downrange  transport  from  the  sixth  floor  varies  from  6  to  15  ft; 
from  the  ground  floor,  3  to  6  ft. 


*  the  top  fragment  in  Figure  25  travels  a  foot  or  so  farther  than  the 
others . 
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Fragment  Distribution,  pgQ  =  30  psi 

The  disposition  of  the  south  and  north  exterior  wall  fragments  re¬ 
sulting  from  the  impact  of  a  30  psi  blast  is  quite  different  from  what  is 
described  above  for  a  1  psi  blast. 

The  time  scale  of  events  is  compressed  approximately  an  order  of 
magnitude  and  wall  debris  is  carried  completely  out  of  the  building. 

Flow  through  the  building  is  established  early  and  the  interior  parti¬ 
tions  being  much  lighter  than  but  offering  the  same  area  as  the  heavy 
masonry  exterior  walls  are  swept  ahead. 

The  southern  external  walls  fail  structurally  at  approximately  6  ms. 

Each  fragment,  top,  side  and  bottom,  is  immediately  driven  into  a 
boundary  of  the  room,  i.e.,  the  floor,  a  partition  or  the  ceiling.  The 
bottom  fragment,  for  example,  is  flat  on  the  floor  41  ms  after  blast 
arrival  and  after  travelling  only  3  ft  downrange.  The  calculations  with 
DEBRIS  show  the  fragment  bouncing  up  to  reach  within  370  ms  a  height  at 
which  the  center  of  mass  is  approximately  2.5  ft  above  the  floor.  Dur¬ 
ing  the  upward  motion  the  fragment  turns  over  slightly  more  than  two 
times  and  travels  downstream  72  ft.  Such  behavior  is  not  realistic  for 
the  fragment  as  a  whole  because  it  will  almost  certainly  break  up  on 
violent  contact  with  the  concrete  floor.  Whether  or  not  the  individual 
masonry  units  composing  the  fragment  will  tend  to  bounce  that  high  or 
not  is  not  known  at  this  time.  However  at  the  peak  of  the  bounce  the 
center  of  mass  has  a  speed  of  41  fps  and  the  blast  wind  is  still  blowing 
strongly.  Considerable  further  downrange  translation  is  to  be  expected, 
although  the  fragment  has  travelled  almost  to  the  north  exterior  wall 
line  during  the  single  bounce. 

At  p  =30  psi  sound  speed  is  1375  fps  so  that  simulated  time  to 
so 

establish  throughflow  at  the  south  facade  after  failure  of  the  north  ex¬ 
terior  walls  is  55  ms.  Translation  of  south  wall  fragments  in  a  classi¬ 
cal  blast  will  commence  at  61  ms;  and  at  the  north  exterior  wall  such 
movement  will  begin  at  172  ms. 


It  is  probably  more  realistic  to  calculate  the  transport  of  the 
bottom  south  wall  fragment  as  due  to  the  sliding  of  individual  masonry 
units  over  the  concrete  floor;  surprisingly,  when  this  is  done  from  the 
point  of  first  fragment  impact  with  the  floor,  and  attributing  the  speed 
of  the  fragment  to  the  individual  unit,  the  net  horizontal  transport  is 
essentially  the  same,  that  is,  the  individual  units  reach  the  north  wall 
line  at  approximately  340  ms.  Their  motion  is  a  tumbling. 

The  top  south  wall  fragment  may  strike  the  ceiling  and  be  deflected 
downward  into  the  floor;  that  possibility  has  not  been  examined  in  a 
DEBRIS  simulation  as  yet.  In  the  absence  of  a  ceiling  the  simulated 
transport  of  the  top  fragment  takes  it  well  past  the  northern  facade,  as 
shown  in  Figure  32.  From  the  sixth  floor,  the  fragment  hits  the  ground 
(at  time  equal  to  2.75  s)  after  the  end  of  the  positive  phase  of  dynamic 
pressure,  and  travels  less  than  10  ft  after  impact.  The  total  downrange 
transport  in  this  example  is  677  ft.  Surface  parameters  supposedly 

O 

appropriate  to  earth  have  been  used,  e.g.,  k  =  0.20  x  10  and  k  = 

10  &hi 
0.35  x  10  dyne/cm.  The  comparable  fragment  in  the  south  wall  of  the 

first  floor  hits  the  ground  outside  the  north  facade  at  approximately 

1.5  s,  which  is  before  the  dynamic  pressure  positive  phase  ends  at  2.67  s, 

but  again,  t'l is  fragment  travels  over  the  ground  only  a  comparably  short 

distance.  The  first  floor  top  fragment  travels  approximately  300  ft 

through  the  air  before  striking  the  ground  and  comer  to  rest  355  ft  from 

its  origin  (Figure  32).  In  these  cases  of  limited  post-impact  transport 

it  is  probably  not  important  to  know  how  the  fragment  breaks  up  on  impact. 

Trajectories  of  the  top  exterior  wall  fragments  from  the  first  and 
sixth  floors  of  the  north  wall  are  depicted  in  Figure  33. 

The  great  contrast  between  Figures  32  and  33  is  surprising.  The 
north  wall  fragments  depart  later  with  a  slower  speed  than  the  south  wall 
fragments  (Table  5),  yet  the  total  downrange  transport  shown  in  Figure  33 
is  twice  that  in  Figure  32.  The  more  slowly  moving  fragment  is  lofted 
higher  and  carried  farther  than  the  initially  fast  moving  fragment. 

The  explanation  lies  in  the  initial  rotation  assumed  for  these  arching 
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exterior  wall  fragments.  The  piece  of  wall  that  breaks  away  early 
with  a  high  speed  rapidly  rotates  out  of  the  attitude  which  is  optimal 
for  lofting  in  the  blast  wind.  The  small  arrows  in  Figure  33  indicate  the 
attitude  of  the  fragment  in  the  locations  shown.  In  Figure  32  the  angular 
speed  is  13  times  its  value  in  Figure  33.  The  realism  of  the  simulation 
in  this  case  remains  to  be  demonstrated. 

As  suggested  above,  the  behavior  of  the  exterior  shear  walls  is 
speculative  at  this  time.  If,  as  seems  likely,  they  fail  during  the 
initial  inward  motion  and  if  they  subsequently  do  not  arch  during  the 
sharp  ir.crease  in  interior  pressure,  they  will  be  flexed  inward  at 
the  time  classical  flow  is  established  through  the  structure.  Such  an 
initial  posture  will  result  in  translation  of  most  wall  debris  down¬ 
stream  and  outward. 

The  progress  of  debris  movement  in  the  Landis  Hospital  at  1  and 
30  psi  is  described  in  Table  6.  Since  at  Pgo=  30  thermal  irradiation  is 
still  underway  during  events  tabulated,  there  is  a  column  to  show  the 
fraction  of  total  irradiation  received  at  each  time  at  the  30  psi  contour 
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Table 

6.  Timing  of  Blast 

Events , 

Landis 

Hospital  (ms) 

Incident 

Ovei pressure  (psi) 

Fraction  of 
Total 
Thermal 
Irradiation 

Event 

1 

30 

at  30  psi 

Blast 

at  south  facade 

0. 

0. 

.350 

South 

exterior  wall  fails 

62. 

6. 

.351 

Blast 

at  I  =  3* 

20.5 

17. 

.353 

Blast 

at  I  =  11 

45. 

25. 

.355 

Collapse  at  I  *  11 

55. 

35. 

.357 

South 

exterior  frags  on  floor 

41. 

.358 

Blast 

at  north  facade 

65. 

55. 

.360 

North 

exterior  fail 

170. 

62. 

.362 

South 

exterior  frags  @  I  = 

3 

170. 

.382 

South 

exterior  frags  @  I  = 

11 

270. 

.401 

South 

exterior  frags  @  I  = 

14 

337. 

.414 

South 

exterior  top  frag  on 

floor 

1000. 

— 

South 

exterior  bottom  frag 

stop 

— 

*  See  Figure  22. 
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Horizontal  Distanc 


CONCLUSIONS  AND  RECOMMENDATIONS 


Over-the-ground  transport  of  heavy  objects  by  blast  can  be  success¬ 
fully  simulated  by  calculation  provided  different  parameter  values  are  chosen 
for  the  three  different  kinds  of  motion,  i.e.,  bouncing,  tumbling  and 
sliding.  In  bouncing  transport  the  object  strikes  the  ground  from  a  height 
and  is  carried  further  downwind;  by  tumbling  motion  is  meant  a  rolling  move¬ 
ment  from  a  position  initially  on  the  ground.  In  sliding  motion  the  object 
does  not  rotate  or  lose  contact  with  the  ground.  The  first  two  kinds  have 
been  quantitatively  observed  in  experiment  and  simulations  are  in  satisfactory 
agreement  with  observations,  although  there  is  insufficient  data  to  develop 
a  relation  between  the  nature  of  the  ground  surface  and  parameter  values. 

The  simulation  BRACOB  has  been  enlarged  and  adapted  to  apply  to  a 
large  urban  building  with  only  modest  increases  in  running  time  on  the  IBM 
3033  computer,  and  the  system  has  been  organized  into  separately  compilable 
modules  so  that  changes  can  be  introduced  with  a  minimum  of  recompilation. 
Capability  to  treat  walls  arching  in  a  frame  has  been  added.  When  combined 
with  the  debris  transport  simulation  and  applied  to  a  weak  walled  building, 
it  appears  to  provide  the  temporal  and  spatial  resolution  required  by  fire 
researchers. 

Just  to  accurately  describe  the  debris  field  in  a  collapsing 
building  requires  improvements  in  the  simulations,  but  perhaps  more 
importantly  at  this  juncture  attention  needs  to  be  given  to  defining 
the  elements  of  the  debris  environment  of  special  importance  to  the 
incidence  and  development  of  fire.  Wall  debris  is  certainly  one  of  the 
important  elements  but  so^  very  likely,  are  debris  arising  in  furnishings 
and  special  air  flow  patterns  within  the  collapsing  building.  This 
may  be  a  good  point  in  this  research  to  introduce  these  considerations. 
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Both  DEBRIS  and  BRACOB  require  improvement.  It  was  originally 
intended  that  fragments  be  of  any  shape  and  there  is  provision  in  the 
DEBUS  code  for  defining  arbitrary  shapes  but  this  feature  has  never 
been  implemented  and  tested.  This  is  most  urgent  in  satisfactorily 
treating  the  so-called  "side"  wall  fragment,  the  trapezoidal  shape  of 
which  may  be  important.  Also  the  code  does  not  handle  without  error  blast 
impact  from  the  side,  i.e.,  an  angle  of  90  degrees.  Some  debugging  is 
needed  in  existing  code.  But  the  most  pressing  deficiency  is  further  data 
for  adjustment  of  the  simulation  parameters,  particulary  to  describe 
the  interaction  with  concrete  surfaces.  Some  further  attention  should 
also  be  devoted  to  additional  checking  of  the  airborne  phase  of  the 
transport  against  existing  data,  such  as  the  small  scale  Bell  Laborator¬ 
ies  models  at  PRAIRIE  FLAT. 

BRACOB  is  in  need  of  enhancement  with  regard  to  the  removal  of 
interior  partitions  which  should  be  automatic  when  the  pressure  differ¬ 
ential  across  the  partition  reaches  a  certain  preset  value. 

After  the  data  from  MILL  RACE  have  been  examined,  both  DEBRIS  and 
BRACOB  might  be  enhanced  to  take  into  account  new  understanding  on 
reinforced  concrete  wall  debris,  establishment  of  flow  through  a 
collapsing  building  and  the  behavior  of  exterior  shear  walls  at  high 
pressures.  BRACOB  presently  drops  calculation  of  wall  response  after  it 
"fails"  structurally;  this  may  be  premature  from  the  point  of  view  of 
debris  distribution.  High  internal  pressures  within  a  room  may  bend  or 
flex  a  "failed"  wall,  in  which  case  its  disposition  as  debris  may  differ 
from  the  present  treatment. 

Finally,  before  the  effects  of  debris  distribution  on  fire  can  be 
understood,  other  categories  of  urban  buildings  must  be  examined  in  the 
way  this  study  has  examined  the  Landis  Hospital.  For  example,  the 
unreinforced  masonry  building  is  important  in  the  eastern  United  States 
and  its  lack  of  a  supporting  framework  will  make  the  behavior  of  the 
floor  structures  different  from  what  has  been  assumed  in  this  study. 
Reinforced  concrete  wall  panels  must  be  treated. 
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Appendix 


CHARACTER  OF  THE  SIMULATED  AIRBLAST 


Original  location  of  the  fragment  is  specified  by  a  free-field 

overpressure  along  with  weapon  yield.  (Height  of  burst  is  assumed 

zero  in  the  work  to  date.)  To  convert  this  input  value  into  range  from 

ground  zero  (GZ) ,  two  methods  have  been  used.  The  first  makes  use  of 

A-l 

the  approximate  relation, 

R  =  |2A2W/[^2Aps  +  B2  -  ^(2Apg  +  B2)2  -  4A2pg2 ^  J1/3  (A-l) 

in  which  side-on  pressure  p  is  measured  in  psi;  radius  R  from  GZ,  in 

1  O  S  £ 

feet;  A  =  3.152  X  10  and  B  =  7.633  X  10  ;  and  W  is  in  megatons. 

Recently,  this  expression  for  R  has  been  evaluated  only  to  obtain 
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starting  values  in  an  iteration  on  the  formula  : 


P 


s 


10.47 
R1 ' 22 


2.9902 


R 


3.053 


4.166 _ 

1  +  0.6096  R 

s 


1.83  “ 


0.2905 


(A-2) 


where  Rg  =  scaled  range,  i.e.,  if  W  =  weapon  yield  in  megatons,  then 

R  =  R/104W1//3 
s 

1/3 

Unit  of  R  is  kft/MT  .  The  differences  between  the  predictions  of 
s 

these  two  methods  are  small,  as  shown  in  Table  A-l.  Neither  formula  is 

used  outside  the  range  of  overpressure  between  1  and  100  psi.  Generali- 

A-2 

zation  of  the  implicit  formula  above  to  any  height  of  burst  is  possible. 

With  original  range  computed  from  given  incident  overpressure,  DEBRIS 
at  each  subsequent  time  step  converts  distance  travelled  by  the  frag¬ 
ment  since  blast  arrival  into  current  range  from  GZ.  If  the  direction 


fwth 


cosine  to  GZ  of  the  y-axis  from  the  original  position  is  c  R  is  oriei- 

y*  o  6 

nal  range,  and  y^  are  original  coordinates,  and  x  and  y  are  current 
fragment  coordinates,  then  current  fragment  range  from  GZ  is  written: 


2  2  2  2 
R  =  (x  -  xq)  +  (y  -  y  )  +  R^  -  2cos(tt 

[(X  -  x±)2  +  (y  -  yo)2]  55 


a!  =  tan  L[(x  -  XQ)/(y  -  y^)] 


a0  =  cos  (c  ) 

2  y 


a  +  a  )R 
1  2  o 


(A-3) 


Although  of  most  importance  to  debris  translation  is  the  character 
of  the  dynamic  pressure  that  is  responsible  for  quasi-steady  drag  of 
a  free  object,  in  some  cases  momentum  imparted  by  the  diffraction  of 
the  shock  front  around  the  fragment  can  be  significant.  When  a  wall 
is  struck  by  blast,  the  diffraction  episode  is  taken  into  account  by 
computing  the  reflected  pressure,  clearing  time,  and  back  pressure 
caused  by  room-filling,  if  any  exists.  The  wall  response  is  computed 
in  terms  of  displacement,  velocity,  and  acceleration  of  an  equivalent 
mass  representing  the  wall  fixed  in  its  edge  supports.  After  the  wall 
is  fragmented  and  is  no  longer  supported,  the  diffraction  episode  is 
usually  over  and  drag  forces  in  the  blast  become  responsible  for  transla¬ 
tion  of  fragments.  Other  objects  of  interest  become  debris  immediately 
upon  blast  arrival,  e.g.,  furniture,  automobiles,  and  the  human  body. 

In  treating  such  objects,  diffraction  momentum  must  be  separately  cal¬ 
culated  and  entered  into  DEBRIS  as  initial  conditions,  i.e.,  as  velocity 
at  a  certain  time  after  blast  arrival.  Displacement  during  the  diffrac¬ 
tion  episode  is  neglected. 

Shock  diffraction  is  a  complicated  phenomenon,  but  good  estimates 
of  its  translatory  effects  upon  initially  free  objects  can  be  made  with 
the  help  of  simplified  concepts.  Essentially  these  amount  to  estimates 
of  front  and  back  pressures,  the  difference  of  which  is  the  instantaneous 
net  pressure.  Consider  a  cube  resting  on  the  ground  struck  head-on  by 
a  blast  wave  of  free-field  static  pressure,  pg.  The  peak  pressure 
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against  the  front  face,  pr,  then  becomes  according  to  standard  analysis 

p  =  2p  (7P  +  4p  )/(7P  +  p  )  (A-4) 

r  s  o  s  o  s 

in  which  P^  =  ambient  pressure,  14.7  psi  at  sea  level.  Initially  the 
back  pressure  is  zero  but  after  a  time  equal  to  shock  transit  time  across 
the  cube,  back  pressure  begins  to  build  up.  During  this  same  time,  front 
face  pressure  is  clearing.  In  fact  to  a  first  approximation  the  dura¬ 
tions  of  these  processes  on  opposite  faces  can  be  considered  equal  to 
each  other  and  approximated  by  the  standard  formula: 

t  =  3S/c  (A-5) 

c 

where  S  =  clearing  distance,  i.e.,  one-half  the  cube  edge  £,  and  c  is 
either  the  sound  speed  in  the  reflected  zone  in  front  or  in  the  ambient 
atmosphere  behind  the  cube.  Front  and  back  pressure  histories  according 
to  these  very  simple  ideas  are  sketched  in  Figure  A-l.  The  diffraction 
momentum,  shown  shaded  in  the  figure,  becomes  the  initial  momentum  for 
the  object  that  starts  moving  at  time  t  =  t  +  tfc  where  tfc  is  the 
transit  time  across  the  cube,  i.e.,  tfc  =  £/U,  where  U  =  shock  speed. 

As  is  evident  in  the  figure  the  complex  pressure  clearing  and  buildup 
processes  In  front  of  and  behind  the  cube  are  reduced  to  linear  changes. 

In  the  case  of  a  cube  1  ft  on  edge  placed  in  the  path  of  a  30-psi 
blast,  initial  momentum  is  estimated  to  be  187  ft-lb/s,  which,  for  a 
cube  weight  of  65  lbs  indicates  a  starting  speed  of  approximately  2.9  ft/s 
at  a  time  equal  to  1.8  ms.  Generally  starting  speeds  of  this  amount  are 
not  important  compared  to  the  speeds  gained  during  the  drag  phase. 

In  the  few  tests  made  so  far,  the  effect  of  including  shock  dif¬ 
fraction  has  been  to  reduce  or  increase  the  total  distance  travelled 
by  the  initially  free  object,  depending  apparently  on  the  original  over¬ 
pressure.  At  high  overpressure  (e.g.,  50  psi),  the  diffraction  episode 
in  effect  substitutes  a  relatively  Inefficient  means  of  locomotion  for 
drag  pressure  during  a  certain  interval  of  time  immediately  after  shock 
arrival.  At  low  overpressure  drag  pressure  declines  dramatically  and 
the  diffraction  of  the  shock  may  be  a  more  effective  means  of  accelera- 
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diffraction 


episode 


ting  the  object.  For  example,  the  simulated  concrete  cube  located 
originally  at  the  50-psi  contour  is  transported  approximately  6  percent 
less  distance  when  it  is  started  at  1.8  ms  with  a  downrange  speed  of 
4.8  ft/s.  In  contrast  accelerating  at  30  psi  by  diffraction  increases 
distance  traveled  by  approximately  6  percent,  and  at  15  psi,  the 
simulated  transport  is  increased  40  percent  by  diffraction.  These 
results  are  suggestive  only  and  some  of  the  effect  may  in  fact  be  caused 
by  other  differences  in  the  calculations. 

No  cognizance  has  been  taken  of  downward  force  of  the  blast  against 
the  ground  springs  through  the  object  resting  on  the  surface.  Presumably 
were  the  surface  perfect,  the  blast  pressure  would  be  effective  in 
compressing  the  spring  with  or  without  the  presence  of  the  object  and 
if  the  surface  is  slightly  rough,  the  static  air  pressure  would  normally 
be  present  under  the  object. 

To  predict  instantaneous  dynamic  pressure,  q(t),  in  this  research, 

two  different  empirical  formulas  have  been  used,  which  generally  differ 

in  their  predictions  by  less  than  10  percent.  The  first  is  based  on 
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simple  curve  fitting  to  the  parameters  in  the  equation  : 


q(t)  =  Qo 


+  (1  -  d)  e_B) 


(A-6) 


where  Q  is  the  peak  incident  dynamic  pressure,  i.e. 
o 

Qq  =  2.5  Pg/(7PQ  +  Ps) 

t+  =  duration  of  the  positive  dynamic  pressure  phase  at  current 
position 

t  =  decay  time  of  blast  at  current  position 

A  =  6  t/t+ 
u 

and 

b  =  <J>  t/t^ 
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The  parameters  d,  6  and  <p  are  given  as  functions  of  static  over¬ 
pressure  p  in  Table  A-2.  For  the  above  formulation  t+  has  also  been 
s  u 

treated  as  a  numerical  parameter  and  its  values  as  a  function  of  over¬ 
pressure  appear  in  Table  A-2  also.  Before  it  may  be  entered  in  the 
formula  for  q(t)  however,  t*  must  be  scaled  to  the  actual  yield  W  (mega¬ 
tons)  : 


+ 

t 

u 


=  10(t+)  .  .  w1/j 

u  tabular 


(A-7) 


The  unit  of  q  is  the  unit  of  p  and  p  ,  but  time  t  must  be  measured  in 

s  o 

seconds.  From  current  fragment  range,  equation  (A-2)  supplies  a  value  of 
incident  peak  overpressure  with  which  to  enter  Table  A-2.  Blast  decay 
time  at  current  fragment  position  must  be  computed  from  shock  travel 
time;  at  each  time  step,  shock  speed  is  computed  from  peak  pressure 
at  current  position  and  divided  into  distance  travelled  during  the 
time  step.  In  this  way  shock  travel  time  from  the  original  position 
to  each  subsequent  position  of  the  fragment  is  available  so  that 
decay  time  at  that  location  is  simply  the  difference  between  current 
time  and  shock  travel  time.  Thus  all  the  quantities  needed  to  evaluate 
q  from  equation  (A-6)  are  available. 

The  second  formulation  of  q(t),  which  is  of  more  recent  origin^ 

makes  use  of  blast  arrival  time  at  the  current  blast  location  in  order  to 

determine  the  interval  t^  during  which  the  blast  has  decayed  at  the 

current  fragment  location,  R.  Corresponding  to  R,  a  value  of  peak  free 

field  peak  overpressure  is  found  as  above,  then  by  iteration  a  blast 
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arrival  time  t  is  determined  : 
a 


(1  +  .6715  t  +  0.004813  t2) 

p  =  i 484 3  - - - y - 

S  (1  +  1.883  t  +  0.02161  t  )  (0.0135  +  t  ) 

a  a  a 


As  a  starting  value  in  this  iteration  t  is  approximated  by 

a 


p  \  -.906088 

*7 


(A-8) 


(A— 9) 
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Table  A-2.  Dynamic  Pressure  Waveshop  Parameters  as  Functions  of  Static 
of  Static  Pressures. 


If  t  is  the  arrival  time  at  the  original  position  and  t  is  the  current 
problem  time,  then  blast  decay  time  at  R  is: 


cd  '  ‘ao  +  ‘  (A-10) 

In  the  equation,  the  computed  arrival  time,  t  ,  has  been  scaled  to 

3 

equivalent  yield  W  (megatons)  of  the  surface  burst. 

The  second  formulation  also  makes  use  of  an  empirical  expression 
for  t  in  terms  of  scaled  range  from  a  1-KT  surface  burst: 

(0.2455  -  0.0115  R  ) 

' _  _ s 

(1  +  61.43  R*) 

and 

Rg  =  R/10(2W)1/3 


(2.177  R  ) 
s 


(1  +  0.7567  R2)+(6 . 147  R3) 
s  s 


» 


0.05546 V  10(2W) 


It  is  the  ratio  of  t^/t^  that  enters  into  the  formula  for  q(t)  as 
follows : 


q(t)  =  Q'o  (R/R’)n 


(A-13) 


As  before  denotes  the  peak  free-field  dynamic  pressure  behind  the 
shock  front,  but  the  primes  make  reference  to  the  state  at  the  current 
shock  front.  Hence, 


Q 


* 

o 


2.5 


7P 

o 


(A- 14) 


in  which  is  computed  from  equation  (A-8)  above  with  arrival  time  equal 

to  t  +  t  ,  scaled  to  1-KT  surface  burst,  i.e.: 
ao 


t’ 

a 


t  +  t 
_ ao 

1/3 

10(2W)  7 


(A-15) 
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Once  is  found,  the  corresponding  radius  R'  may  be  calculated  from 
Equation  (A-2)  above. 

The  numerical  term  n  is  expressed  in  terms  of  current  blast  front 
radius  R'  scaled  to  1  KT:  i.e., 

n  =  0.7917  +  11.04  (R')  +  (14.37  +  6.291  R')/(l  +  28.41  R'3) 

s  s  s 

and 

_  R'  . 

s  (W1/3)  10 

At  the  present  time  the  second  formulation  outlined  above  is  not  appli¬ 
cable  to  explosions  not  on  the  surface. 

Sample  plots  of  instantaneous  dynamic  pressure  histories  prepared 
according  to  the  two  different  formulations  set  forth  above  are  presented 
in  Figure  A-2.  As  can  be  seen,  the  differences  are  small.  Peak  static 
pressure  is  the  parameter  distinguishing  the  several  plots;  however, 
values  of  corresponding  ranges  are  indicated  in  the  figure. 
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Figure  A-2a.  Comparison  of  Instantaneous  Dynamic  Pressure  According  to  Two  Different  Emperical 
Formulations,  P  =  10  psi,  (ranges  correspond  to  W  =  1MT,  HOB  =  0). 


Figure  A-2c.  Comparison  of  Instantaneous  Dynamic  Pressure  According  to  Two  Different  Emperical 
Formulations,  P  =50  psi,  (ranges  correspond  to  W  =  1  MT,  HOB  =  0). 


'  - 


Figure  A-2d.  Comparison  of  Instantaneous  Dynamic  Pressure  According 
Formulations,  p  =90  psi,  (ranges  correspond  to  W  = 
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